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1.0  INTRODUCTION 


The  SEADYN  computer  program  came  into  existence  in  1974  as  an 
extension  of  efforts  to  model  underwater  electronic  structures  at  the 
Electronic  Systems  Division  of  the  General  Electric  Company.  The  major 
features  of  the  program  were  developed  as  part  of  the  author's  doctoral 
studies  at  Cornell  University  (1).  In  1975  the  Naval  Civil  Engineering 
Laboratory  (NCEL)  began  evaluating  SEADYN  and  the  Chesapeake  Division  of 
the  Naval  Facilities  Engineering  Command  funded  extensions  to  include 
moored  vessel  options.  The  output  of  that  funding  was  the  SEADYN/DSSM 
program  and  associated  documentation  (2>S].  Since  1977  the  Naval  Facilities 
Engineering  Command  through  NCEL  has  provided  continuous  development  and 
experimental  verification  of  the  program.  This  activity  has  led  to 
extensive  modifications  of  the  program  and  the  addition  of  new  capabilities. 
These  include: 

1.  A  major  reorganization  of  the  program  structure  and  the 
adoption  of  a  free-field  input  format. 

2.  The  addition  »  '  line  payout/reel-in  capabilities. 

3.  The  addition  of  the  strumming  model  of  Skop,  Griffin  and  Ramberg  (6). 

4.  The  development  of  the  viscous  relaxation  solution  for  static 
analyses  (7). 

The  addition  of  the  bottom-limited  catenary  element. 

6.  The  addition  of  material  internal  damping  models. 

7.  The  addition  of  the  time  sequenced  static  solutions  option. 

R.  The  development  of  a  plotting  post*proces8or  [8]. 

9.  The  addition  of  a  body  impact  model  (developed  for  U.  S.  Coast 
Guard  R&D  Center). 

The  extent  of  these  developments  has  prompted  the  production  of  a 
new  set  of  documentation.  Three  volumes  have  been  written:  the  present 
volume,  a  user's  manual  [9],  and  a  programmer's  reference  manual  (10). 

The  purpose  of  this  manual  is  to  summarize  assumptions,  equations,  and 
numerical  solution  methods  used  in  SEADYN. 
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Finite  element  techniques  based  on  the  stiffness  method  are  employed 
in  SEADYN .  The  hallmark  of  stiffness  methods  is  the  versatility  of 
structural  form  allowed.  Quite  complicated  structural  arrangements, 
loads  and  boundary  conditions  are  permitted.  Unlike  many  special  purpose 
cable  programs,  SEADYN  makes  no  restriction  on  connection  topology  or 
geometric  form.  Multiple  connected  redundant  systems  and  networks  are 
modeled  as  readily  as  single  line  spans.  Since  the  basic  element  is  a 
straight  line  segment,  SEADYN  can  be  used  to  model  truss  structures  as 
well  as  cable  systems. 

Specific  physical  charactf ristics  which  are  important  in  cab?e  and 
mooring  systems  lead  to  nonlinearities  in  the  equations.  Some  of  them 
a  re : 


Geometric  Nonlinearity  -  The  system  stiffness  depends  on  preloading 
.ind  the  deflections  of  the  system. 

Postion  Dependent  Loading  “  Loads  delivered  to  the  system  depend  on 
the  position  and  orientation  of  the  system  and  they  change  as  the  system 
moves.  This  is  typical  of  fluid  induced  loadings  for  drag  and  inertia 
loads  (added  mass). 

Nonlinear  Loading  -  Loads  which  depend  not  only  on  position,  but 
are  a  nonlinear  function  of  the  system  state  variables,  e.g.,  fluid  drag 
loads  depend  on  the  square  of  the  relative  velocity  and  Reynolds'  number. 

Post  ion  Dependent  Constraints  -  The  system  must  remain  within 
specified  constraints,  e.g.,  surface  and  bottom  limits. 

Nonlinear  Materials  -  The  load/strain  relationship  is  dependent  on 
the  amount  of  strain. 

l^hysical  Changes  in  the  Structure  -  The  structure  itself  may  be 
modified  with  time;  snagging  a  body  or  paying  out  lines  are  in  this 
category . 

Each  of  these  phenomena  are  treated  by  SEADYN.  This  manual  outlines 
how  this  is  done. 
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2.0  BASIC  MODELING  ASSUMPTIONS 


The  approach  taken  in  the  SEADYN  computer  prograai  to  model  cable 
and  mooring  systems  can  be  described  as  a  discrete  element  approach.  It 
can  be  considered  as  a  combination  of  the  finite  element  method  and  the 
lumped  parameter  method  in  which  lines  are  modeled  by  the  finite  element 
method  with  bodies  being  lusiped  at  the  node  points. 

Ill  its  classical  form,  the  finite  element  method  seeks  to  represent 
’  I'oiitinuous  physical  systems  with  a  set  of  discrete  or  finite  eleswnts 

which  are  formulated  by  assuming  the  character  of  the  element  response 
in  terms  of  a  set  of  interpolating  functions.  In  its  usual  form  it  is 
equivalent  to  a  Galerkin  form  of  the  method  of  weighted  residuals  where 
the  weighting  functions  are  defined  individually  on  each  element. 

Viewed  from  another  perspective,  the  finite  element  method  is  a  form  of 
the  Kayleigh-Ritz  method  in  which  the  trial  functions  are  defined  only 
on  individual  subregions  (elements)  of  the  system.  The  basic  equations 
for  cables,  mooring  lines  and  hawsers  are  obtained  using  a  sisiple  finite 
element  in  the  form  of  a  straight  line.  The  element  Is  assumed  to  be 
straight  both  before  and  after  deformation  of  the  system,  but  no  restriction 
is  placed  on  the  amount  of  stretch  and/or  rotation  the  element  is  subjected 
to.  It  is  further  assumed  that  bending  and  torsional  effects  are  negligible. 
In  the  case  of  bending  this  means  that  the  bending  stiffness  of  the 
cable  has  negligible  influence  on  the  global  response  of  the  system. 

Neglect  of  the  torsional  effects  does  not  mean  that  twist  is  unianportant. 

It  simply  means  that  coupling  between  twist  and  extension  is  assumed  to 
have  little  effect  on  the  overall  shape  and  response  of  the  system.  An 
obvious  situation  where  this  is  an  invalid  assumption  is  in  low  tension 
conditions  where  a  twist  instability  may  result  in  kinking  or  hockling. 

Ill  addition  to  the  straight  line  element,  a  catenary  model  is  provided 
tor  bottom  interaction.  This  element  neglects  fluid  loads  and  roughly 
jjqif oxiinates  mass  redistribution  with  the  bottom  interaction. 

The  only  deformable  components  in  the  system  are  assumed  to  be  the 
cable  elements.  The  material  is  assumed  to  be  byperelastic,  i.e.,  non** 
linear,  time- independent  with  loading  and  unloading  curves  coincident. 

The  frequency  domain  allows  proportional  damping  effects  while  the  tran¬ 
sient  dynamic  model  allows  various  forms  of  material  damping. 

Any  component  of  the  system  which  cannot  be  modeled  as  an  individual 
line  element  or  a  set  of  line  elements  is  assumed  to  be  a  rigid  body. 

These  rigid  components  are  assumed  to  be  lumped  at  a  single  point  in  the 
system.  They  may  be  assumed  to  have  only  a  point  effect  or  to  act  as  a 
rigid  connector  for  arbitrarily  placed  lines  (e.g.,  a  ship). 

The  system  may  be  totally  immersed  in  a  fluid,  suspended  between  two 
-  fluid.s  (e.g.,  water  and  air)  or  fluid  effects  may  be  ignored.  The  treat¬ 

ment  of  fluid  effects  makes  the  fundamental  assumption  that  the  fluid  and 
struciute  problems  are  uncoupled.  This  means  that  except  for  specific 
localized  effects  the  overall  fluid  field  characteristics  are  unaltered  by 
the  presence  of  the  structure.  Thus,  such  things  as  flow  alteration  due 
to  structural  movement  and  blockage  effects  are  not  dealt  with.  More 
specifics  on  the  assumptions  and  limitations  of  the  fluid  interaction 
_  with  the  structural  system  are  discussed  in  Section  3.5  and  Reference  1. 

3est  Available  Copy 
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Tht*  Lagrangian  approach  ts  taken  in  describing  the  aotion  of  the 
system.  In  this  approach  all  physical  variables  are  expressed  in  terms 
of  their  values  at  an  intitial  reference  state.  It  is  possible  to 
change  the  reference  state  by  employing  generalized  coordinate  transform¬ 
ations  which  account  for  distortions  and  rotations.  Analytical  procedures 
which  begin  from  a  reference  state  and  never  change  that  reference  are 
called  total  Lagrangian.  Updated  Lagrangian  is  an  obvious  title  for 
methods  which  periodically  move  or  update  the  reference  state.  Either 
procedure  ran  be  used  and  the  results  obtained  should  be  equivalent.  In 
the  developments  which  follow  the  configuration  of  a  system  (or  an 
element)  is„designated  by  the  capital  letter  C  and  a  pp-superscript. 

The  symbol  C  means  the  reference  configuration  while  C  aieans  the 
configuration  at  some  time,  t.  The  definition  of  quantities  like  stress, 
strain  and  displacement  usually  involve  two  configurations.  A  pre-su|tcript 
is  used  to  denote  th^  reference  configuration  for  such  cases.  Thus, 
means  a  quantity  in  C  measured  relative  to  C. 

Details  of  the  finite  element  method  applicable  to  cable  systems 
are  given  in  Reference  1.  Only  brief  summaries  of  the  results  pertinent 
to  the  SEADYN  program  are  given  here. 
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3.0  SUMMARY  OF  GOVERNING  EQUATIONS 
3 . 1  Global  Equation  Forms 


The  general  form  of  the  equations  of  motion  for  an  element  can  be 
written 

rIMI  =  Mf)  -  ‘(g)  =  ‘{R}  3-1 

where 

p(M]  is  the  element  mass  matrix 

^{f}  represents  the  external  nodal  forces  in 
^{g}  represents  the  nodal  reactions  in 
^{R]  is  called  the  force  residual 

An  incremental  form  of  the  motion  equations  can  be  written 

J(M1  {Aq}  =  {Af}  -  J(K^1  {Aq}  -  J(C]  {A^}  3-2 

where 


(Aql  =  "^‘(q)  -  ‘(q) 

and 


3-3 


q(K-1  is  called  the  tangent  stiffness  matrix 
K  i 

^[C]  is  an  incremental  damping  matrix 

In  many  situations  (e.g.,  fluid  loading)  the  force  is  dependent  on 
the  deflection.  In  this  case 

(M)  =  ijM  it  t  (Aq)  3-4 

3  J(q) 

=  (if)  +  J(Kjl  (iq) 
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The  incremental  motion  equations  are  then  written 


J(M1  {Aq’l  +  J(C1  {Aq}  ♦  ^|K^|  {Aq}  =  {A?}  3-5 

where 

r(Kt)  = 

Equation  (3-2)  or  (3-6)  can  be  used  to  model  small  displacement 
response  about  a  steady  deformed  configuration,  or  it  can  be  used  in 
nonlinear  dynamics  by  recalculating  the  stiffness  matrix  at  each  step. 

It  should  be  noted  that  both  equations  neglect  the  position  dependent 
effects  in  the  mass  matrix.  The  incremental  load  rotation  matrix,  iKj.), 
is  nonsymmetric  and  causes  some  problems  in  applying  Equation  (3-6). 

Its  effect  when  small  increments  are  used  is  felt  to  be  minimal  and  is 
ignored  in  SEAOYN.  The  incremental  damping  matrix  siay  be  difficult  to 
obtain  in  the  more  general  situations.  A  simplified  treataaent  is  discussed 
in  Section  3.3. 

An  alternative  form  of  the  incremental  equations  is  obtained  from 
Equation  (3-1)  by  expanding  only  the  internal  loads  in  a  Taylor  series 
and  neglecting  higher  order  terms.  The  result  is 

♦  Rlltfl  (Aq)  « 

This  equation  is  linearized  by  approximating  ^nd 

^ith  their  values  at  t*^At  while  remaining  in  the  ori|^^£tion  defined  by 
^C.  Any  damping  effects  are  assumed  to  included  in  '^{f}. 

The  contributions  from  each  of  the  elements  in  the  system  (cables, 
lumped  bodies,  and  rigid  bodies)  can  be  combined  in  a  very  siaq>le  and 
direct  manner  once  they  have  been  generated  and  transformed  to  the 
global  coordinate  system.  This  is  done  element  by  element  by  accumlating 
the  element  contributions  in  the  appropriate  position  of  the  global 
arrays.  An  ordering  of  the  degrees  of  freedom  is  implied  in  this  procedure. 
The  order  assumed  in  SEADYN  is  simply  the  three  displacement  components 
(x,y,z)  stored  in  the  order  of  the  node  number.  Thus,  the  global  nodal 
displacement  vector  assumes  the  x  component  of  node  number  one  is  first, 
the  y  component  of  node  2  is  fourth,  etc.  By  requiring  slave  (movement 
defined  in  terms  of  another  node)  nodes  to  be  numbered  after  nodes  %diich 
have  active  degrees  of  freedom,  the  solution  bookkeeping  is  greatly 
simplified . 

The  assembled  global  equations  have  essentially  the  same  form  as 
the  element  motion  equations.  The  main  distinction  is  that  the  order  of 
the  equations  is  increased  to  include  all  of  the  active  degrees  of 
freedom  in  the  system.  Noting  this,  the  total  nonlinear  equations  of 
motion  can  be  written 
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JiMJ  J{q}  ^  *{f}  -  ^{g}  =  *'{RJ 

The  two  increoiental  forns  are 

rINI  (A<i)  *  J  |C|  lAql  ■>  JlK^I  (Aq)  =  (Ml  3-9 

"“hi  ♦  rIKtI  IAq|  =  "“Ul  -  ‘t*l  3-10 

it  should  be  emphasized  that  these  represent  the  assembled  equations 
for  the  system  and  that  it  is  assumed  that  the  constraints  iaq>lied  by 
the  boundary  conditions  and  slave/mastcr  conditions  are  accounted  for. 

The  dynamic  equations  reduce  to  the  static  equations  when  the  time 


dependent 

terms  are 

dropped.  Thus,  the  nonlinear 

static  equation  is 

*-{R} 

=  0 

3-11 

and  the  incremental 

static  equations  are 

j^IKtI 

{Aq}  = 

lAfl 

3-12 

r(Kt] 

{Aq}  = 

"“ifl  -  Mg) 

3-13 

In  the  siatic  case  the  parameter  t  is  used  to  signify  a  load  step  rather 
than  a  time  step.  The  static  equations  presume  a  stable  physical  system 
has  been  described  by  imposing  adequate  constraints  on  the  system.  If 
this  is  not  so,  the  stiffness  matrices  are  singular  and  there  is  not  a 
unique  configuration  of  the  system  which  will  satisfy  the  equation.  In 
ease  of  nonlinear  systems  (particularly  those  with  surface  ships  and 
mooring  buoys)  the  static  global  equations  may  be  ill-conditioned.  This 
means  they  are  nearly  singular  and  numerical  errors  in  the  solution  pro¬ 
cedure  may  lead  to  apparent  singularities.  This  will  be  given  more 
attention  in  Section  4.2. 

The  classical  approach  for  analyzing  wave  induced  motions  of  plat¬ 
forms  and  vessels  is  to  transform  the  incremental  equations  into  the 
frequency  domain  and  use  linear  superposition  techniques.  The  linearised 
small  displacement  equation  form  used  is  that  of  Equation  (3-9).  The 
transforming  assumption  is  that 

{M}  -  3-14 

Assuming  quasi-linearity  of  Equation  (3-9) •  the  steady  state  reapottse 
has  the  form 
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Subst i tut ilig  (3*  14)  (3-lS)  into  (3”9)  yields 


|M)  +  iu) 


^  JlK 


!Q} 


{F} 
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Solution  of  this  set  of  complex  simultaneous  linear  equations 
allows  the  computation  of  the  response  amplitudes  and  phase  angles  for 
all  degrees  of  freedom  in  the  system.  The  magnitude  of  the  response  in 
each  degree  of  freedom  is  given  by 


3-17 


where  Qi  is  the  complex  conjugate  of  the  i  component  of  {Q} .  The 
ph.se  angle  between  the  incident  loading  (wave)  and  the  response  is 
gi  en  by 


=  tan 


1  ’m 
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))hase  angle  of  zero  corresponds  to  the  response  being  in  phase  with 
the  incident  loading. 


?. . 2  F.quations  for  Line  Eleaients 


Two  line  elo'iiont  models  are  available.  The  first  is  the  one-diroen- 
sional  simplex  element.  This  is  a  straight  clement  using  two  nodes. 

The  second  element  also  uses  two  nodes  in  its  definition  and  is  in  the 
torm  of  a  catenary. 

3.2.1  The  One  Dimeusioiial  Simplex  Element 

A  finite  element  which  has  the  form  of  a  line  (i.e.,  one-dimensiona 
and  ii.sos  only  the  field  parameters  at  the  two  ends  of  the  line  in  the 
interpolating  function  is  referred  to  as  a  one-dimensional  simplex 
eleme*nt  (11 j.  When  the  element  is  a  straight  line  in  3D  space  and  the 
field  parameters  are  the  nodal  (end  point)  displacements  the  element  is 
called  a  truss  element  in  structiiral  terminology. 

Consider  a  single  straight  element  which  is  defined  by  the  position 
of  two  nodes  (one  at  each  end).  Select  a  local  coordinate  system  with 
the  X  axis  extending  from  the  first  node  to  the  second.  The  other  two 
axes  may  be  chosen  arbitrarily  under  the  restriction  that  they  form  a 
right  lianoed  carte.si3n  reference  frame.  When  the  element  is  in  its 
unloaded  state  it  has  a  length  I..  Assume  the  material  constitutive 
relation  has  the  form 


o 


t 


F, 


o 


t 

c 

o 
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where 


^'S  is  the  2nd  Piola-Kirchhoff  stress  in 
o 


is  the  Green's  strain  in 
o 


E  is  a  nonlinear  material  modulus  which  may  be  a  function  of  strain. 


The  incremental  form  of  t^^s  constitutive  relation  can  be  written  for  a 
small  strain  increment,  Ao 


t+Atg  ^  *^S  -f 

o  o  0  T  o 
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.  tE..  is  the  tangent  modulus  evaluated  at  t. 
where  o  T  *  C. 


Making  the  finite  element  assumption  that  the  displacement  at  any 
position  along  the  element  ?s  a  linear  function  of  the  displacements  of 
the  nodes  one  can  write 


{>*] 
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whe  re 


{u}  represents  the  components  of  the  displacement  from  "c 
{q}  represents  the  components  of  the  nodal  displacements 
^[Nj  is  called  a  shape  function  matrix 
is  the  identity  matrix  of  order  3 

kinematic  relations 


(Green 


The  symbolic  expression  for  the  j^^tge  displacement  kinemati* 
n's  Strain)  for  a  movement  from  C  to  C  can  be  written  [l] 


^{c}  -  JlD]J{uJ 
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SuhsHtutiou  from  EquuLiun  (3-21J  yields 


‘(.I  =  JiDI  ,|N|  JUI  =  J  181  J(,l 


The  mass  matrix  for  the  straight  element  can  be  written  in  two  forms: 


Consistent  Mass  Matrix 


|Ml  = 


Luinped  Mass  Matrix 


K  Ra  Rt 

iMi  =  -e_4-JL 
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R  R  R  R 

where  p  is  the  element  material^density  in  C,  and  C,  and  A  is  the 

element  cross-sectional  area  in  C.  With  the  assumption  of  conversion 

of  mass,  the  element  mass  matrix  does  not  change  with  deformation.  The 

pre-sub  and  superscripts  are  used  in  Equation  (3-1)  since  this  is  not 

true  of  fluid  added  mass. 

The  consistent  mass  matrix  is  obtained  from  the  kinetic  energy  and 
Equation  (3-21).  The  lumped  form  can  be  obtained  by  the  intuitive  pro¬ 
cess  of  lumping  half  of  the  element  mass  at  each  node  or  by  summing  all 
the  terms  on  each  row  of  the  consistent  mass  matrix  and  assigning  the 
sum  to  the  diagonal  position. 

t 

The  external  forces  may  be  due  to  point  or  distributed  loads. 

Point  loads  appear  as  specific  entries  in  the  global  equations.  Distributed 
loads  are  usually  from  gravity  effects  and/or  fluid  loading.  Fluid 
loading  effects  are  discussed  in  Section  3.5.  The  general  form  of  the 
gravity  loading  is 


"{f}  = 


J'  Rjj^jT 
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where 


p 

{T}  represents  the  components  ol  the  element  specific  weight 
(in  fluid)  relative  to  the  local  coordinate  system 


^(Nj^  is  the  transpose  of  ^(Nj 
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Substitution  from  Equation  (3-21)  into  Equation  (3-26)  and  noting 
the  orientation  of  the  element  with  respect  to  the  direction  of  gravity 
leads  to  the  conclusion  that  these  forces  are  equivalent  to  placing  one 
half  of  the  element  weight  acting  in  the  gravity  direction  at  each  node. 
It  should  be  noted  that  Equation  (3-26)  assumes  mass  is  conserved. 

The  internal  forces  of  Equation  (3-1)  can  be  written  (1} 


3-27 


where  {A)  is  a  un^^t  vector  in  the  direction  of  the  deformed  element. 
The  stress  term,  j^S,  can  be  written 

t-q  - 

where  is  the  element  load  in  the  deformed  state.  This  allows  the 
force  residual  to  be  written 


‘'{Ri  =  ‘'{f}  -  ‘'P 
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The  stiffness  matrix  can  be  obtaiu|>d  from  consideration  of  the  second 
variation  of  the  strain  energy  in  '"C.  The  result  is  [1] 


j^(Kt1 
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where 


Ik 


Ik^l 


o  T 


o^-T  * 


1  0  0 
0  0  0 
0  0  0 
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1^21 


I^gI 


1  0  0 
0  1  0 
0  0  1 


/ 


r”2 


r”i 


3.2.2  The  Bottom  Limited  Catenary  Element 

Quite  often  a  mooring  system  will  employ  a  line  which  roust  interact 
with  the  bottom  as  the  moored  body  moves.  The  major  feature  of  such  a 
line  i.s  that  significant  lengths  of  line  are  lifted  or  laid  down  on  the 
bottom  when  relatively  small  movements  are  induced  at  the  upper  end  of 
the  line.  This  is  particularly  true  of  shallow  water  moors.  As  an  aid 
ill  modeling  such  line.s  SEADYN  provides  a  bottom-limited  catenary  element. 

The  element  uses  the  classical  catenary  relations  with  a  modifica¬ 
tion  to  account  for  length  changes  due  to  line  stretch.  The  nature  of 
the  catenary  equations  is  such  that  it  is  not  possible  to  solve  explicitly 
for  stiffness  terms.  A  perturbation  procedure  presented  by  Peyrot  and 
Goulois  (12)  is  used  to  circumvent  this  problem.  The  SEADYN  implementa¬ 
tion  follows  the  developments  of  Reference  12  with  two  exceptions.  The 
first  is  that  the  tangential  stiffness  matrix  produced  by  the  perturbation 
procedure  wa.s  recognized  to  be  deficient  since  it  did  not  recognize  the 
geometric  stiffness  resisting  out  of  plane  motion.  This  Jeficiency  was 
removed  by  adding  a  factor  to  the  diagonal  of  the  tanger.t  stiffness 
matrix  at  the  degree  of  freedom  corresponding  the  out-of-plane  displacement 
of  each  node.  The  factor  added  is  the  horizontal  component  of  the 
catenary  tension  divided  by  the  horizontal  distance  between  the  nodes. 

The  second  modification  consists  of  the  imposition  of  a  bottom 
limit  on  the  catenary.  The  limit  is  assumed  to  be  a  horizontal  plane 
passing  through  the  lowest  node  on  the  catenary.  The  approach  taken  to 
enforce  this  limit  is  an  iterative  process  similar  to  the  one  proposed 
by  Peyrot  113]. 
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The  actual  equations  used  ace  those  of  Reference  12  and  will  not  be 
repeated  here.  The  procedure  relys  on  an  algoritho  which  computes  the 
nodal  forces  on  a  sagged  line  when  the  unstretched  length,  nodal  positions 
and  the  distributed  loading  given.  The  tangent  stiffness  natrix  is 
obtained  from  the  changes  in  nodal  loads  produced  by  perturbing  the 
nodal  positions.  Only  cable  weight  loading  is  treated.  Fluid  drag 
loads  are  ignored  on  this  element. 


The  bottom-limited  catenary  element  uses  a  local  coordinate  system 
which  has  its  origin  at  the  bottom  node.  The  local  x  axis  is  horizontal 
in  the  vertical  plane  containing  the  two  element  nodes.  The  positive  x 
direction  is  toward  the  upper  node.  The  local  y  axis  is  vertical  upward. 
The  z  axis  is  chosen  to  form  a  right-handed  cartesian  system.  Element 
stiffness  and  contributions  to  the  residual  are  first  computed  in  the 
local  coordinate  system  and  then  transformed  to  the  global  system.  The 
transformation  procedures  are  those  discussed  in  Section  3.9. 

3 . 3  Material  Models,  Elasticity 

The  line  material  is  assumed  to  be  nonlinear  clastic  in  form.  Two 
functional  forms  are  assumed  which  relate  the  line  tension  to  the  exten- 
.sional  strain.  The  first  is  a  tabular  form  which  represents  a  sequence 
of  linear  segments  describing  the  relation.  The  second  is  a  two  parameter 
curve  fit  form.  It  has  been  found  effective  in  modeling  cable 
constructions  Il3].  The  form  is 


where  i  and  b  are  curve  fitting  parameters.  Note  that  this  relation  is 
assi.mei  to  be  between  the  total  tension  and  the  total  uniaxial  Green' s 
strain.  The  relation  between  the  customary  uniaxial  engineering  strain 
and  uni.ixial  Green's  strain  is 


-  (1  +  i  t  .  )  f.  3-32 

v>  2  engineering  engineering 

The  diife  ence  between  these  two  strain  measures  is  minor  up  to  strains 
oi  iM  an  engineering  strain  of  20%  the  Green's  strain  is  only  22%. 


The  tr.aterial  const i t.utive  relations  are  used  in  two  situations  in 
.SKf\P''N,  The  first  occasion  arises  in  the  case  where  an  initial  equili- 
brii.i.i  (.'onf igurati on  is  known  (or  guessed).  In  this  case  it  is  desired 
to  f.nd  the  unstrained  length  and  the  strain  with  the  load  and  position 
of  thr'  nodes  given. 


In  the  other  case,  the  unstretched  liriglh  is  known,  the  nodal  posi- 
tio  i:  are^given  and  it  is  desired  to  cofni. -e  the  strains  and  the  tangent 
modulus,  ®A.  This  situation  occurs  time  a  new  estimate  of  the 

state  of  ?he  system  is  made. 
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In  either  case,  it  is  necessary  to  use  the  expression  for  uniaxial 
Green's  strain 


=  1 

2 

1  '  ^ 

0  2 

This  form  is  used  rather  than  accumulating  strain  increments  to  avoid 
numerical  round-off  errors.  The  for  each  element  is  either  gi^en  or 
computed  for  each  element  at  the  beginning  of  the  analysis.  The  is 
recomputed  at  each  stage  of  the  solution  process  where  the  strain  is 
needed.  Th£  computation  of  these  two  lengths  can  be  a  source  of  numerical 
noise  when  ~  “L  and  low  precision  is  carried  in  the  computer.  Two 
situations  aggravate  the  problem:  (1)  the  material  is  very  stiff,  i.e., 
large  (2)  the  loading  is  small. 

3.4  Damping  Models 

Damping  enters  the  SEADYN  analyses  in  various  ways.  First,  there 
is  a  certain  amount  of  damping  inherent  or  intentionally  included  in  the 
numerical  algorithms  for  calculating  transient  res,  onse.  Such  algorithmic 
damping  is  treated  in  Section  5.2.  Another  form  of  numerical  damping  is 
involved  in  the  nonlinear  static  solutions.  This  is  an  artifice  to  aid 
in  obtaining  convergent  static  solutions  to  highly  nonlinear  problems. 
Discussion  of  numerical  damping  c.an  be  found  in  Section  5.1.  A  third 
form  of  damping  is  due  to  the  dissipative  effects  of  the  fluid  interacting 
with  the  structure  as  discussed  in  .Section  3. .5. 

I'his  section  is  concerned  with  material  or  physical  damping  models 
for  the  dynamic  response  of  the  structure.  Material  damping  is  considered 
to  be  an  inherent  char.acteristic  of  the  line  material  and  is  distributed 
throughout  the  line.  It  is  appropriately  treated  in  the  material  consti¬ 
tutive  relation.  The  result  is  that  the  material  takes  a  viscoelastic 
form.  An  alternative  formulation  presumes  not  enough  information  is 
avail.ible  to  define  all  features  and  sources  of  dissipation  and  the 
behavior  is  approx iin.ated  through  an  estimated  damping  matrix.  The  form 
assumed  is 

J;1C1  =  a  ^(M]  +  P^(K|  3-34 

The  treatment  of  material  damping  tollows  a  component  approach. 

Two  components  are  used  as  building  blocks:  the  elastic  component  which 
relates  forces  to  displacements  and  the  dashpot  which  relates  forces  to 
velocities . 

These  components  may  be  used  in  two  arrangements  in  SEADYN.  These 
are  referred  to  as  the  Kelvin  and  the  NOAA-Reid  models  (251.  They  are 
repr€!seiited  in  Figure  3-1. 
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a)  Kelvin  Model 


b)  NOAA'Reid 


These  relations  can  be  written  in  terms  of  strains  using  the  substitutions 


K6  -  EAe 

n6  =  -^6  =  CAe 

Then 

Kelvin  Mode 1 

T  =  EA  e  +  C  A  e  3-37 

o 

NOAA  -  Reid  Model 

CA  CA- 

T  =  EA  e  +  ^  (EA  +  EA.)  i  *  f  3-38 

o  EA^  o  1  EAj 

Equation  (3-37)  has  been  implemented  in  SEADYN  by  recognizing  the 
EA  term  as  the  one  discussed  in  Section  3.3.  The  CA  term  is  taken  as 
an  additional  input  parameter. 

The  more  complicated  form  of  equation  (3-38)  presents  more  problems 
.since  it  involves  a  tension  rate  term  as  well  a.s  a  straiu  rate  term.  As 
with  the  Kelvin  model,  the  EA  term  is  identified  as  the  basic  mate  ial 
elasticity  relation  and  the  cR,  and  EAj  parameters  are  provided  as  addi¬ 
tional  input.  The  tension  rate  is  approximated  by  a  backward  difference 
in  time.  As  a  result  it  lags  the  other  parameters.  This  approximation 
i.s  considered  to  be  a  reasonable  compromise. 

These  material  damping  models  have  only  been  implemented  in  SEADYN 
on  the  direct  iterative  (D1 )  method  for  transient  dynamics  (see  Section  5.2). 
This  solution  works  directly  with  the  global  equations  in  the  form  of 
Equation  (3*8)  and  no  global  stiffness  or  dumping  inatrix  is  used.  The 
1)1  method  makes  the  use  of  proportional  damping  a  little  less  direct. 

In  The  case  of  proportional  damping  in  transient  ■  namics  using  the  DI 
method  it  is  assumed  that  for  each  element 


/  0,  2  ^  „  t,,  o.  ^  t  • 

(o(  m  1,  +  8  t  A)  c 

^  o  s  o 


whe  re 


o  s 


is  the  material  secant  modulus. 


3-39 


The  frequency  domain  solutions  treat  damping  by  forming  an  incremental 
damping  matrix  corresponding  to  the  form  of  Equation  (3-16).  The  material 
damping  forms  have  not  been  implemented,  but  the  proportional  damping  in 
the  form  of  Equation  (3-34)  has.  (ii  addition,  the  dissipative  terms 
from  drag  loading  have  been  approximated  using  the  approach  presented  in 
Appendix  1. 
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Neither  inateii<il  uor  pruporLiunal  damping  is  used  in  any  transient 
dynamic  option  besides  the  D1  method.  No  damping  is  used  in  the  natural 
frequency  calculations. 

3 .  .I  Flu  id  Loads 

The  primary  assiuaption  regarding  the  effect  of  fluid  immersion  on 
the  cable  system  as  stated  in  Section  2  is  that  the  fluid  and  structural 
problems  are  uncoupled.  >'he  independence  principle  [l4]  is  assumed  in 
the  treatment  of  fluid  loading  on  the  cables,  In  brief,  the  independence 
principle  asserts  that  the  fluid  loading  can  be  treated  as  resulting 
from  two  separate  flows:  one  normal  and  one  tangential  to  the  cable. 

The  fluid  induced  loads  can  be  separated  into  a  part  involving  the 
relative  velocity  between  the  cable  and  the  fluid  and  a  part  involving 
the  relative  acceleration.  The  velocity  related  terms  can  be  written 


{w}  =  Wjj  In}  + 


3-40 


who  re 


=  1/2  D 


w.j,  ^  1/2  Pj.  C.J,  D  V,^ 

In). 


('  C 
N’  T 


are  unit  vectors  in  the  directions  of  the  normal 
and  tangential  components  of  the  relative  velocity 

are  normal  and  tangent i.il  drag  coefficients 

is  the  drag  diaoieter  of  tl»e  cable 

is  the  <lensity  of  the  fluid 

are  normal  <u»d  laugenl »al  relative  velocities 
The  tliiid  loading  vector  for  an  clement  is  then 

^f}  =  f  {w}  dl. 
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linear  variation  of 
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The  acceleration  related  portion  of  the  fluid  loading  can  be  separated 
into  a  part  due  to  flow  field  acceleration  and  a  part  due  to  structural 
acceleration.  The  flow  field  acceleration  part  is  neglected  in  the 
cable  loading  and  the  structural  acceleration  part  Is  treated  as  an 
added  mass.  The  added  mass  matrix  for  a  straight  element  can  be  written 


^^added^ 


r  «  Rt 
2 


0 
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where 


is  an  added  mass  coefficient. 

It  should  be  noted  that  there  is  no  added  mass  tangential  to  the  cable 
which  makes  this  a  position  dependent  term. 

In  some  ot  the  solution  procedures  economies  can  be  achieved  if  the 
mass  matrix  is  diagonal.  A  lumped  form  of  the  combined  cable  and  fluid 
added  mass  matrix  can  be  written: 


(M)  = 


Cu  P,  A  1. 


2  2 


where 


iMj  t 


P  = 


pt 


3-46 


3-47 


This  form  still  presents  some  difficulties  which  will  be  dealt  with  in 
the  discussion  of  the  solution  procedures. 

Fluid  loads  on  submerged  lumped  bodies  are  estimated  using  an 
approach  similar  to  that  used  for  the  cable  clement.  Two  forms  for 
lumped  bodies  are  considered:  spherical  and  an  ond-faired  cylinder. 
The  drag  loading  on  a  sphere  is  given  by 

^{f  .  }  =  1/2  p^  Cp,  n  *{A}  3- 

sphere’  'f  D 


where 

D  is  the  diameter  of  the  sphere 
Cjj  is  the  drag  coefficient 
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V  is  lh«‘  r«!lativr  v«*locity  t  li<‘  fluid  and  the  point 

whore  tho  si>hero  is  locat^.-d. 

Mm  is  a  unit  vector  in  the  direction  of  relative  velocity. 
The  added  mass  for  a  sphere  is 


•”added^  "" 


C  p  n  IV 

_rj — i - 1 1  I 

6  ‘  ‘3' 


The  end  faired  cylinder  loading  and  added  mass  i.s  assumed  to  have 
the  same  form  as  that  for  a  cable  element  except  that  the  total  effect 
is  placed  at  a  .single  point  rather  than  being  distributed  between  two 
nodes . 

There  are  specific  fluid  effects  pectiltar  to  surface  buoys  and  ships 
which  are  discussed  in  Sections  3.7  and  3.8. 

!1.6  [.limped  llodiej 

I'wo  forms  of  lumped  bodies  are  considered.  The  simplest  form 
tre.it.s  the  body  as  a  single  point  with  three  di spl acemeiil  degrees  of 
freedom.  The  point  is  assumed  to  have  mass  hut  no  rotational  inertia. 

No  el.isticity  eflocts  are  attributed  to  the  body,  hut  it  m.ay  be  a  means 
of  inducing  fluid  loads  into  the  system  (see  Section  T-.S).  If  the  mass 
of  the  body  is  m,  then  the  mas.s  matrix  is  .simply 


I  M  I  I  4  ~  I" 

lumped 


Tins  body  affects  only  one  node  in  t  h«*  system. 

The  .second  form  is  that  of  a  rigid  body  with  spatial  dimensions. 

The  nia.s.s  is  still  assumed  to  bt.-  concentrated  at  a  single  point  but  that 
point  has  six  degrees  of  freedom.  I’he  ma.s.s  matrix  assumed  for  this  case 
lias  specific  forms  only  in  the  case  of  mooring  buoys  and  surface  ships. 
The  program  generates  the  mass  matrix  for  mooring  buoys  and  assumes  the 
mass  iii.'itrix  is  defined  by  input  from  the  ship's  motion  file  for  .surface 
i.liips.  See  Sections  3.7  and  3.8  for  details.  In  either  case  the  mass 
matrix  is  a  b  x  (>  matrix. 


3.7  Surface  Ships^  and  I’^at forms 

The  rigid  body  element  is  used  to  model  ships  and  platforms  (the 
term  ship  will  be  used  to  mean  either  one).  A  single  node  point  is  used 
to  define  the  position  of  the  .ship.  Since  the  node  must  express  the 
angular  as  well  as  spatial  position  of  the  ship  it  is  reipiired  to  have 


ix  degrees  of  freedom.  (The  SEADYN  progr.am  uses  two  consecutive  nodes 
of  t,hi'«?<>  degrees  of  freedom  each  to  define  a  ship.)  .Mtaclunents  of 
mooring  lines  and/or  working  lines  are  handled  through  the  slavc/master 
(  raiisformation. 


In  static  analyses  the  steady  state  effects  of  winds  and  surface 
currents  acting  on  the  ship  are  treated  as  lateral  and  longitudinal 
forces  and  a  yaw  moment  which  are  assumed  to  act  at  the  ship's  reference 
point.  The  values  of  these  forces  depend  on  the  flow  velocities  and  the 
angle  between  the  ship’s  heading  and  the  flow  direction.  Empirical  load 
tables  giving  load  coefficients  versus  heading  or  analytical  load  functions 
may  be  used.  The  empirical  approach  is  given  in  NAVFAC  DM-26  (15]  and 
both  approaches  are  summarized  in  tie  appendices  of  the  user's  manual. 

The  dynamic  equations  for  the  response  of  surface  ships  to  waves 
are  usually  given  in  an  incremental  linearized  form.  These  equations 
have  the  form 

t^^S  ^  ^AS^  ^*^8^  *  t^^S^  *  t^^S^ 

where 

{Ugl  represents  the  six  components  of  ship's  motion 
(surge,  sway,  heave,  roll,  pitch,  yaw) 

^[M-J  is  the  ships  mass  matrix  including  rotational  inertial  terms 

L  O 

is  the  added  mass  due  to  fluid  acceleration  effects 

t,  As> 

Mc„]  is  an  equivalent  linerarized  damping  matrix 

t.  o 

^[KgJ  is  the  ship's  hydrostatic  restoring  matrix 

{fg}  are  the  point  equivalent  forces  representing  the  wave 
induced  exciting  forces 

In  order  to  obtain  this  linearization  it  is  usually  assumed  that 
the  ship  is  driven  by  a  simple  harmonic  wave.  With  this  assumption,  the 
forces,  added  mass  and  damping  are  frequency  and  heading  dependent.  In 
addition,  the  linearization  of  the  roll  damping  term  makes  it  dependent 
on  the  magnitude  of  the  roll  angle.  Equations  of  the  form  of  (3-51)  can 
be  obtained  for  slender  bodies  using  strip  theory  (16],  A  more  general 
theory  is  required  for  other  forms  (17).  The  SEADYN  program  assumes  the 
value.s  lor  these  coefficients  are  provided  through  a  data  file.  The 
format  of  that  data  file  is  described  in  Appendix  A  of  the  User's 
Manual  and  Section  5.3  of  the  Programmer's  Reference  Manual. 

The  element  equation  represented  by  Equation  (3-51)  can  be  manipulated 
as  any  other  element  equation  and  combined  with  the  global  equations  of 
the  system.  No  new  concepts  are  involved  in  these  operations. 

J.8  Mooring  Buoys 

Ship's  moors  often  involve  surface  buoys  which  support  the  aM>oring 
line  and  are  connected  to  the  ship  through  a  hawser.  This  type  of  buoy 
usually  remains  on  the  surface  where  it  is  subjected  to  the  effects  of 
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winds,  currents  and  waves.  The  general  fors  of  the  buoy  aotion  equations 
linearized  to  represent  snail  excursions  fron  a  static  reference  state 
can  be  written: 


t 

t 


t 

t 


(Kggl  {Ugg} 
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where  the  various  terms  follow  the  previously  established  pattern.  The 
forcing  term  represented  by  {fen}  deals  only  with  wave  excitation. 

Static  load  effects  on  a  mooring  buoy  follow  the  sasie  form  as  described 
previously  in  Section  3.6. 

In  order  to  avoid  a  grt:at  deal  of  complexity,  it  is  assumed  that 
the  buoy  is  spherical  in  shape  and  that  a  local  cartesian  coordinate 
system  is  selected  which  is  vertical  in  the  z  direction  and  has  the 
incident  wave  traveling  in  the  *x  direction.  No  loss  of  generality  is 
incurred  with  this  choice  of  coordinate  system  on  a  spherical  buoy  since 
the  character  of  the  coefficients  do  not  depend  on  the  orientation  or 
attitude  of  the  buoy.  The  problem  is  further  siiaplifled  if  it  is  assumed 
that  the  buoy  is  homogeneous  with  the  center  of  gravity  at  the  geometric 
center  of  the  buoy  and  that  the  geometric  center  is  located  at  the  water 
line  in  the  reference  state. 


Attachments  of  hawser  and  mooring  lines  to  the  buoy  can  be  readily 
handled  if  their  positions  relative  to  the  local  coordinate  system  are 
known.  The  rigid  link  transformation  described  in  Section  3.11  is  used 
for  this  purpose.  The  positions  of  the  attachments  can  be  found  from 
the  static  solution. 

Since  it  is  assumed  that  the  equations  use  as  the  reference  con¬ 
figuration,  the  configuration  notation  will  be  dropped  for  this  discussion. 


Transformation  from  the  local  to  the  global  system  for  assembly  of 
the  buoy  equations  with  the  rest  of  the  moor  system  equations  is  a 
straightforward  process  which  follows  the  method  outlined  previously. 

For  this  reason  only  the  coefficients  in  the  local  coordinate  systm 
will  be  given  here.  It  should  be  kept  in  mind  that  the  following  equations 
represent  the  incremental  motion  equations  for  a  surface  buoy  in  the 
local  coordinate  system  just  described. 


Given  that  the  buoy  has  a  mass  designated  by  m  and  a  mass 
inertia,  J  ,  the  buoy  portion  of  the  laass  matrix  is 


m 


IHjI  = 


m  I. 


nt  of 
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where  I.  is  the  identity  matrix  of  order  3.  The  assumption  of  a  homoge¬ 
neous  sphere  should  be  recalled  at  this  point.  If  the  attacluwnts 
contribute  significant  bmss,  their  effects  can  be  treated  by  including 
additional  lumped  masses  at  those  nodes. 
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The  added  mass  has  the  form 


A  0 

XX 


"xO 

0  0 


0  0  0 


(SYM) 


‘“sir  ^  ‘“b'  "■  ‘“a^ 

The  wave  damping  matrix  has  a  form  similar  to  the  added  mass  matrix. 
Specific  values  for  the  added  mass  and  damping  coefficients  for  a  sphere 
were  given  by  Patton  [20].  llis  value.s  were  obtained  by  curve  fitting 
the  analytical  results  presented  by  Kim  (18].  The  nondimensional  values 
obtained  were 


M  =  1.089  0.052  a 

X 


for  0  <  a  <  0.79 


]/(-0. 0.118  +  0.959  a  )  for  0.74  <  a  <  .3.4 


M  =  1.85 

Y. 


1 .02  a 


'-0.256 


for  0  <  a  <0.1 


0. 1  <  a  <3.4 


.V  -  0 

X 


for  0  <  a  <0.1 


=  -O.069  t  0.715  a  for  0.1  <  a  <  1.37 


for  1.37  <  a  <  3.4 


0.126  +  1.7  a 


f  o  r  0  <  a  <  0 . 4 


1  .  Ifie 


-0.83a 


f o  r  0 . 4  <  a  <3-4 


whore  a'  is  the  nondimensional  frequency  given  by 


2/Ta  _  a  u) 
^  ~  g 


22 


and 

a  =  radius  of  sphere 
A  =  wavelength 

w  ~  circular  frequency  of  wave  container 

These  equations  for  the  added  mass  and  damping  coefficients  for 
mooring  buoys  proved  to  be  cumbersome  for  programming.  A  polynomial 
curve-fitting  of  the  original  curves  given  by  Kim  [18]  was  used  in  the 
program.  The  polynomial  coefficients  are  given  below. 

FUNCTION  a^  a^  a^.  a^ 

M  1.0620-0.4090  5.3299-10.0143  7.6387-2.9089  0.5509-0.0414 

X 

M  1.7945  1.3362  -10.9227  18.0521  -13.8920  5.5725  -1.1249  0.0902 

z 

N  0.  -1.5252  7.2144  -8.6447  5.0535  -1.6218  0.2750  -0.0192 

X 

N  0.  4.3747  -9.8378  10.7232  -6.6657  2.3676  -0.4439  0.0339 

z 
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and  V  =  the  kinematic  viscosity  of  the  f^iiid^  Since  the  kinematic 
viscosity  of  water  is  of  the  order  of  10^  ft  /sec,  the  rotational  added 
mass  terms  will  be  small  compared  to  the  buoy  inertia  terms.  Damping 
due  to  rotational  motion  is  very  small  and  will  be  neglected.  Thus, 


When  the  center  of  pressure  does  not  coincide  with  the  center  of 
gravity  of  the  buoy,  a  coupling  between  lateral  and  rotational  motion 
exists.  These  terms  can  be  written: 


xe 

A  (z  -  z  ) 

XX  eg  cp 
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y4' 

-A  (z  -  z  ) 

yy  ‘-'g  cp^ 

'xO 

C  (z  -  z  ) 

XX  eg  cp 
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-y4( 

-C  (z  -  z  ) 

yy  eg  cp 

With  the  origin  of  the  local  coordinate  system  at  the  geometric 
center  (also  center  of  gravity),  z,  is  zero.  The  center  of  pre.ssure 
for  a  tialf  submerged  sphere  is  obtained  from 


e  z  r]  dS 
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The  damping  terms  presented  above  do  not  represent  the  effects  of 
v:  Kcou.s  drag.  The  viscous  terms  involve  the  square  of  the  relative 
velocity  between  the  buoy  and  the  fluid  and  are  therefore  nonlinear. 
The  VISCOUS  effects  are  generally  of  less  import. ince  than  the  wave 
damping.  Obviously,  this  is  not  the  case  for  the  rotational  movement 
since  those  terms  are  zero  for  v.ave  damping.  in  a  free-floating  buoy 
the  vi.scous  rotational  terms  would  play  an  important  part,  but  in  a 
mooring  system  where  the  hawser  and  mooring  leg  restrain  the  buoy  the 
rotation  is  limited.  Therefore,  all  of  the  viscous  terms  will  be 
neglected  rather  than  attempting  to  l.inearize  them. 


The  only  nonzero  hydrostatic  restoring  force  on  a  hal f-submerged 
spherical  buoy  acts  in  the  heave  direction.  Its  value  for  small  displace¬ 
ments  is 


*^22 


na^  pg 
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When  the  buoy  provides  a  connection  between  a  mooring  line  and  a  hawser 
it  develops  an  additional  stiffness  (resistance  to  motion)  due  to  the 
tensile  force  being  transmitted  across  it.  This  is  analogous  to  the 
geometric  stiffness  term,  seen  in  the  cable  element  stiffness 

matrix  in  Section  3.2.1.  Thxs  geometric  stiffness  effect  is  automatically 
taken  into  account  by  the  rigid  link  transformation  of  the  elements 
representing  the  attached  lines. 

The  right-hand  side  of  Equation  (3-S2)  represents  the  forces  due  to 
surface  waves.  Assuming  the  wave  is  harmonic  in  form,  Kim  [19]  shows 
that  the  wave  induced  forces  can  be  written 
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where 


pa^  U) 


^  (  a'(z+ix)  „ 

pg'^g  e  n 


^  f  a'(z+ix) 
pg-^S  e 


ds  ~  pg-'g  cos  (a'x)  ds 

(xn  -  zn  )  ds 
'z  'x 


For  the  half  submerged  sphere 

A^  S  _  ,  -  0 

^00  “  00  ^2  ”  ° 


The  wave  pressure  component  of  the  heave  exciting  force,  1^,  is 
closely  approximated  by  the  pressure  at  the  water  surface  distributed 
over  the  cross  section  at  the  water  surface.  A  plot  of  this  function 
versus  the  nondimens ignal^freguency  js  given  in  Figure  3-1. 

The  values  for  M  ,  w  and  N;  for  a  half  submerged  sphere  were 
by  Kim  (Ref  19).  A  polynomiaf  curve^flt  of  those  functions  plus  the 
curve  for  are  summarized  in  the  following  table. 


SUMMARY  OF  WAVE  EXCITING  FORCE  COEFFICIENTS 


FUNCTION 

"o 

"i 

"2 

"3 

*4 

“5 

^6 

X 

0. 

1.7586 

-8.2171 

12.0253 

-8.2882 

3.0081 

-0.5576 

0.0416 

z 

1.7868 

1.0552  - 

10.5792 

17.2729 

-13.2596 

5.3272 

-1.0781 

0.0866 

N® 

X 

1.0833 

0.0833 

- 1 . 4496 

0.7753 

-0.0705 

-0.0314 

0.0055 

0. 

n" 

z 

0. 

4.2382 

-8.5367 

8.2743 

-4.6849 

1.5633 

-0,2823 

0.0211 

1 j/na^pg 

1.0 

-0.0004 

-0.1218 

-0.0026 

0.0069 

-0.0007 

0. 

0. 

where  the  function  form 

is : 

F  = 

N 

L 

a  (a' )^ 

3' 
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Thus  all  of  the  terms  necessary  for  treating  the  small  displacement 
behavior  of  a  restrained,  half-submerged  spherical  buoy  are  available. 
The  assumptions  employed  appear  to  be  reasonable  and  should  lead  to  a 
good  approximation  of  the  buoy  effects  in  a  deep  sea  moor.  Although 
a  spherical  buoy  has  been  assumed,  a  comparison  of  the  curves  presented 
by  Kim  (Ref  18)  for  a  sphere  and  those  presented  by  Garrison  (Ref  17) 
for  a  ha  If -submerged  cylindrical  buoy  with  an  aspect  ratio  of  1.0  shows 
that  the  added  mass  and  damping  coefficients  are  quite  similar.  There¬ 
fore,  it  is  reasonable  to  expect  the  sphere  equations  to  give  at  least 
an  order-of-magnitude  approximation  of  a  cylindrical  buoy. 

The  small  displacement  assumption  deserves  some  further  comment. 

For  wavelengths  of  the  order  of  the  buoy  diameter  one  would  expect  the 
buoy  motion  to  be  small.  However,  as  the  wavelength  increases,  the 
buoy  motions  increase.  Since  most  of  the  wave  energy  is  expected  in 
the  longer  wavelengtits ,  the  buoy  could  be  expected  to  see  large  motions 
which  would  cause  these  equations  to  be  inaccurate.  When  no  ship  is 
in  the  system,  this  inaccuracy  could  be  serious.  Fortunately,  the  ship 
motion  becomes  a  significant  effect  in  the  longer  i./avelengtbs  and  the 
buoy  motion  is  dominated  by  the  ship  movement  as  transmitted  through 
the  hawser  and  reacted  by  the  mooring  line.  In  this  case  the  contri¬ 
butions  from  the  buoy  itself  (though  in  error)  would  generally  be 
insignificant.  For  the  shorter  wavelengths,  the  ship  appears  nearly 
fixed,  and  the  exciting  forces  due  to  wave  action  on  the  buoy  become 
important.  It  is  fortuitous  that  this  is  the  range  in  which  the  buoy 
equations  are  roost  accurate. 

3 . 9  Coordinate  Systems  and  Transforma tions 

Equations  for  the  cable  element,  as  well  as  ships  and  mooring 
buoys,  are  most  readily  developed  in  a  local  or  intrinsic  coordinte 
system  which  is  considered  to  move  with  the  clement.  The  development 
of  3  global  set  of  equations  which  represents  the  behavior  of  the 
assembled  system  of  elements  requires  a  single  global  coordinate  system. 
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Transformations  between  the  local  coordinate  system  and  the  global 
system  in  a  reference  configuration  are  accomplished  by  the  usual 
rotation  of  coordinates.  The  general  form  for  the  components  of  a 
vector  at  a  point  is 


jITl 


global 


3-78 


Columns  of  the  transformation  matrix  are  the  components  of  a  unit  vector 
in  the  direction  of  the  local  coordinate  axis  expressed  in  the  global 
system.  Since  only  right-handed  cartesian  coordinates  are  considered, 
the  inverse  transformation  is  obtained  with  the  transpose  of  ij[T) .  All 
contributions  from  the  external  and  internal  loads,  etc.,  must  be 
transformed  to  the  global  system  before  they  are  combined  with  other 
components  in  the  system. 

It  should  be  noted  that  these  coordinate  transformations  apply  only 
to  quantities  expressed  relative  to  a  specific  reference  configuration. 
As  long  as  the  reference  configuration  remains  fixed  the  individual 
coordinate  transformations  Remain  ^nchanged  regardless  of  how  much 
deformation  occurs  between  c  and  C. 


Coordinate  transformations  involving  nodal  point  displacements  for 
an  element  can  be  written 


local 
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7he  form  for  nodal  forces  is  similar.  When  one  transforms  the  incremental 
motion  equations  from  the  local  to  the  global  system  the  equations  take 
the  form 


^IMI 

^^^^global  ^  r‘^1  ^^^^global 

+  =  ^^^^global 
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Thus,  it  is  seen  that  coordinate  rotations  do  not  alter  the  form  of 
the  equations.  Unless  specific  emphasis  is  required,  no  further  distinc¬ 
tion  between  the  local  and  global  equation  forms  will  be  made.  It  is 
assumed  that  the  equations  are  written  in  a  homogeneous  system,  i.e.,  all 
displacements,  forces,  stiffnesses,  etc.,  are  in  the  same  coordinate  system. 

3.10  Restraint  Transformations 


A  generalization  of  the  coordinate  system  transformation  is  useful 
in  modeling  the  effects  of  displacement  restraint.^ .  Such  restraints 
represent  boundary  conditions  where  the  value  of  the  displacement  is 
specified.  It  may  be  zero  or  some  finite  quantity.  SEADYN  imposes 
these  restraints  in  static  analyses  after  the  global  stiffness  matrix 
has  been  assembled.  The  process  amounts  to  an  imposition  of  a  linear 
relation  of  the  form: 


3-81 


Here  it  is  assumed  the  (<!]  represents  known  or  specified  displacement 
components. 


Consider  a  partitioned  form  of  the  static  stifiiiess  equations  (the 
indices  are  dropped  for  simplicity) 
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Applying  the  linear  transformation  represented  by  (3-81)  the  following 
results 


'“.J  '“.I  = 
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The  form  of  this  equation  suggests  a  procedure  for  imposing  restraints 
on  the  global  equations.  First  assume  the  d-portion  has  only  one  component, 
then  take  the  column  of  the  stiffness  matrix  that  corresponds  to  that 
degree  of  freedom,  and  subtract  it  times  d  from  the  force  vector.  Next, 
s<^t  .ill  entries  in  that  row  and  column  of  the  stiffness  matrix  to  zero. 
Finally,  replace  the  entry  in  the  force  vector  for  that  row  to  the  value 
of  d  and  set  the  corresponding  diagonal  element  of  the  stiffness  matrix 
to  unity.  This  process  is  repeated  for  each  degree  of  freedom  where  a 
restraint  is  to  be  imposed.  This  process  is  valid  for  any  value  of  d 
(including  zero).  When  d  represents  a  large  movement  it  will  be  necessary 
to  utilize  some  nonli.iear  solution  algorithm  to  get  an  equilibrium  solu¬ 
tion.  These  methods  ..re  discussed  in  Section  5.1. 
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3 ^  ^  Slave/Mastcr  Transformations 

The  generalized  rigid  bodies  used  to  model  ships  and  mooring  buoys 
require  some  special  manipulations  to  connect  them  to  the  rest  of  the 
system.  Their  motions  are  assumed  to  be  described  by  the  six  degrees  of 
freedom  of  a  single  point.  The  att.achmsnts ,  however,  will  not  connect 
to  the  rigid  body  at  that  node.  Since  the  body  is  assumed  to  be  rigid, 
it  i.s  possible  to  express  the  motion  of  any  point  on  the  body  in  terms 
of  the  motion  of  a  single  point  and  the  relative  positions  on  the  body. 
The  node  used  to  model  the  body  is  called  a  master  node.  Any  other 
point  on  the  body  is  called  a  slave  node.  Given  the  six  components  of 
motion  at  the  master  node,  the  translation  components  of  a  slave  node 
can  be  written: 
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master 

or 

‘‘^^slave  ~  ^*^^naster 

•Wiere  Ax,  Ay,  Az  are  components  of  the  distance  between  the  two  points 
measured  from  the  master  to  the  slave,  i.e.,  Aa:  = 

The  matrix  (Tgul  t)e  viewed  as  a  generalized  form  of  the  coordinate 
transformation  represented  by  Equation  (3-7ft)  and  the  transformation 
procedures  of  the  previous  sections  can  be  employed. 

It  should  be  noted  that  the  slave/niaster  transformation  involves  an 
alteration  of  the  number  of  degrees  of  freedom.  When  an  end  of  a  cable 
element  connects  to  a  slave  node,  the  application  of  Equation  (3“8A)  in 
the  transformation  indicated  by  Equation  O'SC)  results  in  a  stiffness 
matrix  (etc.)  which  is  9  x  9  instead  of  6x6.  It  should  further  be 
noted  that  the  slave/master  transformation  form  assumes  small  displace¬ 
ments  and  is  therefore  only  applicable  to  the  incremental  equations. 
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where  hl'T)  is  a  rotation  gatrix^of  the  form  (3-78)  which  represents  the 
total  angle  changes  from  u  to  C  at  the  master  node. 
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4.0  MODELING  CADLK  AND  MOORING  SYSTEMS 
4 . I  General  Approach 

The  SEADYN  program  utilizes  a  line  element  in  such  a  way  that  it 
can  represent  truss  structures  as  well  a;;  cables  and  mooring  lines.  All 
ot  these  structures  have  a  common  feature  in  that  the  basic  element  is  a 
one-force  member.  Modeling  a  cable  span  with  multiple  line  segments  or 
representing  a  truss  structure  involves  essentially  the  same  steps.  Of 
course »  a  segmented  model  of  a  curved  span  represents  some  special 
approximations  not  needed  for  trusses  and  the  geometric  nonlinearities 
must  be  treated  in  cable  spans  while  that  may  not  be  essential  in  trusses. 
But  once  the  approximations  have  been  n»ade,  the  building  block,  the 
element,  is  one  that  can  be  used  for  either  type  of  structure. 

One  should  be  aware  that  whenever  only  two  line  elements  meet  at  a 
node  and  the  other  ends  of  those  elements  are  not  totally  fixed,  the 
geometric  nonlinearity  must  be  treated.  V/ithout  preloading  and  appropriate 
external  loads,  such  a  system  is  a  mechanisn..  The  stiffness  matrix  for  a 
mechanism  is  singular. 

Although  truss  structures  can  be  defined  that  have  sufficient 
rigidity  and  structural  stability  that  small  defurmatiua  assumptions  can 
give  reasonable  solutions,  small  deformation  theory  cannot  be  used  on 
cable  spans  unless  large  preloads  are  supported  and  the  stiffening 
effects  of  those  preloads  are  included  in  the  equilibrium  equations. 

In  general,  the  modeling  process  for  cable  systems  involves  the 
selection  of  appropriate  subdivision  of  the  spans.  Since  straight  lines 
are  the  pror'.ominaut  elements  in  SEADYN,  it  is  necessary  to  make  the  sub¬ 
division  fine  enough  to  capture  the  curvatures  of  the  initial  and  final 
st.ates  with  atcoptable  error.  Here  curvature  requires  more  elements. 

Other  modeling  assumptions  relate  to  material  variations  and  lumped  body 
approximations.  Nodes  must  be  located  where  materials  change  (only  one 
material  per  element)  and  where  bodies  are  to  be  lumped.  Nodes  al.so 
must  be  placed  where  limit  conditions  or  other  restraint/constraint  con¬ 
ditions  may  be  .sp<>cifiod.  Liberal  use  of  nodes  and  elements  has  various 
economic  impacts  and  good  modeling  practice  leads  to  rational  compromises 
between  economics  and  accuracy. 

Proper  modeling  procedures  requires  some  special  understanding  of 
the  loading  environment  and  the  possible  boundary  conditions.  It  is 
necessary  to  separate  loading  conditions  into  categories.  Three  categories 
dealt  witl)  in  SEADYN  are 

DEAD  -  Loads  inherent  with  the  structure  (weight,  buoyancy) 

which  have  virtually  no  temporal  variation.  They  must 
be  supported  along  with  all  and  any  other  loads. 

LIVE  -  Loads  which  are  temporary  in  nature  but  applied  slowly 

enough  to  be  assumed  static.  Current  and  wind  loads 
can  he  pl.iced  in  thi;;  category.  Various  operational  o; 
working  loads  are  of  this  type. 
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Transient  -  Loads  which  excite  or  produce  acceleration  effects 
leadinfi;  to  movement  in  a  time  scale  short  enough  to 
preclude  static  approximations.  Generally  this  will 
mean  oscillatory  responses  in  which  mass  related  forces 
are  not  negligible. 

SKADYN  assumes  all  three  of  these  conditions  will  bo  encountered 
and  since  geometric  nonlinearities  are  involved,  the  behavior  in  one 
condition  is  likely  to  be  dependent  on  the  state  induced  by  other  conditions. 
A  typical  situation  is  one  in  which  dynamic  responses  arc  highly  sensitive 
to  the  static  state  produced  by  .1  combination  of  giavity/bouiidary  loads 
and  current  loading  with  negligible  time  variations  but  significant 
space  variations.  It  is  not  appropriate  to  attempt  dynamic  analy.sis  in 
this  case  until  the  static  solution  is  found.  The  structure  of  SKAUYN 
presumes  a  staged  or  sequenced  analysis  process  which  applies  new  loading 
conditions  to  the  state  obtained  in  the  proceeding  stage.  Full  account 
is  given  for  the  geometric  and  property  changes  in  this  staging.  A 
typical  analysis  sequence  would  be  to  compute  the  equilibrium  state  for 
DEAD  loads,  apply  the  LIVE  loads  and  finally  induce  the  transient  loads 
on  the  combined  st.itic  state. 

The  primary  variables  in  the  computations  are  the  nodal  displace- 
nioiits  and  nodal  velocitie.s.  Secondary  information  such  as  strains  and 
tensions  arc  comtuted  from  the  primary  data.  Interpretation  of  results 
should  always  he  vempered  with  an  understaiiding  of  the  approximations 
and  methods  used.  Since  the  straight  elements  pre.suiiie  linear  functions 
for  displ. a  cements  between  the  nodes  only  a  constant  value  for  element 
strain  can  be  obtained.  The  conscqiicnce  of  this  is  a  single  value  for 
the  element  tension, 

Care  is  needed  in  specifying  initial  data  for  highly  flexible 
.structure's  with  low  prelo.ad.  Such  problems  siifiei'  from  a  tendency  to 
i-.'idically  change  shape  when  .small  variations  are  made  in  element  lengtlis, 
initial  node  positions  or  element  tensions.  One  often  encounters  great 
difficulty  in  obtaining  numerically  .st.iblc  initial  configurations  for 
dead  loading.  More  is  said  of  this  in  later  sections  and  the  various 
references . 


.  2  S  lu 

1  Mooring  Systems 

The 

surtaci*  !■ 

basic  components  required  in  modeling  mooring 
ihips  are: 

.systems  for 

a . 

surface  ship 

b. 

mooring  iine.s  (usually  submerged) 

c . 

mooring  buoys 

d. 

hawser  (usually  in  air  and  subjected  to  wind 

lo.'ids) 

e  . 

floats,  sinkers  and  anchors 
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Each  of  these  comi)<>nent s  have  l.een  dealt  with  to  some  extent  in  the 
previous  discussion.  Unfortunately  it  is  not  possible  to  develop  a 
fully  K^neral  nonlinear  analysis  of  ship's  mooring  dynamics.  The  primary 
reason  for  this  is  the  highly  complex  nature  of  the  interaction  between 
the  sea  and  surface  bodies  such  as  ships  and  mooring  buoys.  The  equations 
presented  for  these  bodies  are  linearized  equations  which  address  only 
the  response  to  harmonic,  long-crested  waves. 

The  theoretical  approach  used  in  dealing  with  mooring  systems 
follows  through  a  series  of  approximatious .  The  first  step  is  to  obtain 
a  description  of  the  mooring  system  in  the  quie.scent  state  where  only 
gravity  loads  are  involved.  This  is  called  a  dead  load  analysis.  There 
are  some  pitfalls  in  this  step  which  are  related  to  the  geometric  nonlin¬ 
earity  of  the  problem.  One  does  not  usually  know,  a  priori,  what  the 
dead  load  configuration  of  the  system  is  and  the  terms  in  the  static 
equations  are  configuration  dependent.  This  is  the  so-called  initial 
configuration  problem  and  it  is  dealt  with  in  some  detail  in  References  1 
and  7.  Section  4.8  also  di.scusses  the  problem. 

The  active  loads  on  the  mooring  system  are  assumed  to  be  winds, 
surface  and  subs^irface  currents  and  surface  waves.  The  effect  of  surface 
waves  is  primarily  a  dynamic  phenomenon  while  the  static  or  steady-state 
effects  of  winds  and  currents  predominate.  Therefore,  it  is  assumed 
tha(  Use  active  loads  on  the  system  can  be  separated  into  a  static 
eltect,  which  is  primarily  due  to  wind.s  and  current.s,  and  a  dynamic 
effect  from  surface  wave.s. 

The  next  stop  in  the  analysi.s  is  then  a  static  analysis  to  obtain 
the  response  to  winds  and  currents.  Thi.s  is  referred  to  as  a  live  lo'*d 
analysis  and  it  may  include  point  loads  representative  of  imposed  work 
load.s.  Nonlinearilies  in  the  sysiein  also  play  an  important  part  in  a 
live  load  analy.sis.  Tlie  geoinolric  nonlinearity  is  still  present.  A 
significant  new  nonlinearity  loiue.s  from  the  noncouservat  ive ,  position 
dependent  loads.  The  flow  induced  loads  are  .strongly  dependent  on  the 
orient.ition  of  the  various  elemouts  and  the  Ig.-i-Is  change  direction  and 
magnitude  as  the  system  changes  shape  and  position. 

M.iterial  noni  ineurit  ies  also  play  an  import.ant  part  at  this  stage. 

The  most  pronounced  effect  ccines  from  the  fact  that  the  lines  cannot 
.support  compiessi ve  loads.  Should  the  imposed  loads  cause  any  of  the 
legs  to  go  slack  (usually  in  taut  moors  with  neutrally  buoyant  lines), 
the  material  stiftiiess  goes  to  zero  and  the  group  of  elements  involved 
with  the  log  are  part  of  an  unstable  structure  and  the  stiffness  matrix 
becomes  singular. 

An  iaiportaiit  feature  in  slack  moors  (negatively  buoyant  legs)  is 
the]'-  ability  to  resist  loads  by  changing  shape  and  by  lifting  line  or 
laying  it  down  at  the  bottom.  Modeling  this  interaction  with  the  bottom 
using  reasonably  long  cable  elements  produces  some  approximations  that 
must  be  ki'pt  in  mind  while  interpreting  the  results.  An  alternative  is 
to  use  the  bottom  limited  catenary  eU'ment,  This  element  also  introduces 
appruximat ions  that  must  be  considered. 
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In  certain  situations  it  is  not  possible  (or  feasible)  to  obtain  a 
dead  load  configuration  and  then  proceed  to  the  live  load  analysis. 

Direct  solution  for  the  combined  dead  and  live  loads  is  usually  required 
in  these  cases.  One  example  is  the  solution  for  a  single  point  moor. 

In  this  case  the  quiescent  state  is  of  little  value  and  solution  for  the 

combined  effects  is  usually  quite  easy  to  obtain  if  one  has  a  good  esti¬ 

mate  of  the  total  horizontal  load  to  be  supported  by  the  moor.  Another 
situation  where  this  procedure  may  be  required  is  in  dealing  with  taut 
moors  where  legs  go  slack.  In  this  case  the  slack  legs  could  be  ignored 

(not  included  in  the  model)  and  the  dead  load  included  in  the  live  load 

analysis . 

The  analysis  of  wave  induced  dynamics  begins  with  the  static  configu¬ 
ration  developed  by  the  wind  and  currents.  Of  necessity  this  is  a  fre¬ 
quency  domain  analysis.  The  equations  presented  for  the  ship  and  mooring 
buoy  dynamics  were  obtained  by  assuming  the  excitation  was  from  a  harmonic 
wave.  The  linearization  process  renders  the  equations  frequency  dependent 
and  limited  to  small  motion  amplitudes. 

The  harmonic  loading  input  assumes  a  reference  point  which  is  the 
defining  node  for  the  ship.  Phase  angles  on  the  loading  are  induced  in 
the  loading  which  are  dependent  on  the  harmonic  wave  length  and  the  dis¬ 
tance  of  the  load  point  from  the  reference  point.  Response  phase  angles 
represent  shifts  from  loading  applied  at  the  reference  point. 

The  dynamic  tension  response  in  each  element  is  obtained  by  adding 
the  displacement  increments  to  the  nodal  positions  recalculating  the 
element  tensions  from  the  constitutive  relations,  using  the  new  lengths, 
and  subtracting  the  tensions  in  the  static  reference  state. 

Wave  loading  produces  mean  and  very  low  frequency  forces  which  may 
have  significant  magnitudes  on  moored  vessels.  These  depend  on  the 
responses  and  orientation  of  the  vessel  relative  to  the  wave.  These 
wave  induced  drift  forces  lead  to  slowly  varying  loads  on  the  system 
which  lead  to  shifts  in  the  "static”  reference  position.  The  treatment 
of  these  forces  requires  at  least  two  passes  through  the  frequency 
domain  solution.  After  the  first,  the  magnitudes  of  the  drift  forces 
are  estimated.  Since  these  are  predominately  static  loads,  the  static 
analysis  should  be  repeated  with  the  additional  loading.  If  this  results 
in  significant  movement,  the  frequency  domain  solution  should  be  repeated 
on  the  new  reference  state. 

4.3  Statics 


Various  load  incrementing  options  are  provided  with  the  static 
analyses.  In  general,  the  gravity  load  will  be  increased  in  increments 
in  a  dead  load  analysis.  Fluid  induced  loads  are  increased  in  incresients 
during  live  load  analyses.  Both  analyses  allow  point  loads  which  are 
incremented  from  zero  to  the  specified  value,  held  constant,  or  reduced 
from  the  value  to  zero.  It  is  usually  assumed  that  gravity  loads  are 
held  fixed  during  a  live  load  analysis,  but  an  option  is  provided  which 
allows  gravity  load  to  be  incremented  during  a  live  load  analysis. 
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Wlien  ships  are  used  in  a  static  analysis  the  ship  is  constrained  to 
remain  on  the  surface.  This  normally  means  the  ship  is  fixed  in  heave, 
roll  and  pitch.  It  is  possible  to  input  stiffness  terms  to  allow  small 
heave  and  angular  motions  out  of  the  surface  plane.  No  checks  are  made 
on  the  magnitude  of  those  responses,  so  the  user  is  warned  to  evaluate 
the  validity  of  his  results  when  these  stiffnesses  are  used. 

^ ^  Dynamics,  Time  Vs.  Frequency  Domain 

SEADYN  provides  computation  options  for  dynamics  in  the  time  and 
frequency  formats.  Transient  dynamic  analysis  in  the  time  domain  presumes 
all  of  the  nonlinearities  are  active.  It  is  further  assumed  that  only 
the  line  elements  provide  equations  for  the  elastic  strain  energy. 

Kinetic  energy  terms  (mass  related  terms)  are  obtained  from  the  lines 
and  lumped  bodies.  Dissipative  terms  come  from  material  daaq)iDg  and  the 
drag  loading.  No  equations  are  available  within  SEADYN  for  the  rigid 
bodies.  These  typically  require  major  computations  using  fluid/solid 
interaction  equations  which  are  beyond  the  scope  of  SEADYN. 

A  time  domain  solution  is  called  for  whenever  the  problem  imposes 
motions/ loads  which  result  in  large  movements,  load  variations  producing 
nonlinear  material  effects,  compression  of  flexible  lines  (slack), 
and/or  modifications  of  limit  conditions.  The  time  domain  form  is  the 
only  one  capable  of  dealing  with  those  nonlinearities.  The  restrictions 
mentioned  above  on  rigid  bodies  eliminate  a  very  important  area  of 
computation  from  transient  analysis  in  SEADYN.  On  the  positive  side, 

SEADYN  provides  a  significant  capability  for  analysis  of  difficult 
dynamic  problems:  anchor  deployment;  array  placement,  adjustment  and 
response;  towing  of  lines  and  bodies  in  irregular  paths;  mooring  leg 
response  when  attachment  point  motion  is  known,  etc. 

The  frequency  domain  solution  is  a  quasi-linear  approach  aimed 
directly  at  the  moored  ship/platform/buoy  problem  discussed  in  Section  4.2. 
It  is  applicable  when  the  assumptions  of  DEAD,  LIVE  and  wave  loading  are 
appropriate.  It  is  significant  in  that  it  provides  the  ability  to  treat 
coupled  responses  between  the  mooring  lines  and  the  moored  rigid  bodies. 
Highly  restricted  models  for  spherical  mooring  buoys  are  built  into 
SEADYN,  but  the  body  equations  and  loading  functions  must  be  provided 
externally  for  ship/platform  components. 

The  time  domain  and  frequency  domain  forms  can  be  used  in  tandem  to 
investigate  mooring  leg  dynamics.  The  approach  would  be  to  follow  the 
combined  mooring  solution  through  the  frequency  response  t(»  define  the 
motion  at  the  mooring  line  attachment.  This  motion  is  then  converted  to 
a  time  sequenced  description  which  is  used  to  drive  the  top  of  a  stati* 
cally  preloaded  mooring  leg.  This  produces  a  more  realistic  definition 
of  the  vessel  motion  and  then  allows  assessment  of  the  mooring  leg  non- 
linearities. 
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4  5  Natural  Frequencies  and  Mode  Shapes 

SEADYN  includes  the  option  for  cosaputing  natural  frequencies  and 
mode  shapes  simply  for  informational  purposes.  It  provides  data  to  no 
other  option  within  SEADYN  when  the  global  structure  is  investigated. 

The  computation  of  strumming  induced  drag  amplification  uses  this  option 
to  get  strum  string  mode  shapes  on  subsets  of  the  structure. 

The  computations  for  the  global  structure  take  the  lumped  (diagonal) 
mass  matrix  without  correction  for  the  lack  of  tangential  added  mass 
along  with  the  tangent  stiffness  matrix  for  currently  defined  state. 

The  tangential  added  mass  correction  is  not  taken  since  it  leads  to  a 
non-diagonal  mass  matrix,  and  that  requires  additional  storage.  Co,..puter 
storage  is  a  serious  issue  on  this  option  since  the  Jacobi  method  is 
used.  The  Jacobi  method  works  entirely  in  core  and  requires  the  storage 
of  a  full  nxn  matrix. 

4 . 6  Surface  and  Bottom  Constraints 

The  foregoing  discussion  introduces  the  problem  of  constraining  the 
model  of  the  system  to  responses  which  lie  between  the  natural  boundaries 
imposed  by  the  water  surface  and  the  bottom.  One  does  not  expect  buoys 
or  lines  to  rise  out  of  the  water,  or  lines  and  anchors  to  go  below  the 
bottom.  The  SEADYN  program  assumes  the  surface  and  bottom  are  flat  and 
parallel.  Checks  ace  made  at  each  step  of  static  and  time  domain  analyses 
to  see  if  nodes  of  the  system  are  within  the  imposed  limits.  To  avoid 
the  costly  operation  of  checking  all  nodes  in  the  system  at  every  step, 
checking  is  only  done  on  points  where  special  limit  conditions  are 
spocifieJ.  When  one  of  the  critical  nodes  is  within  a  certain  tolerance 
distance  of  the  surface  or  bottom,  it  is  constrained.  For  buoys  or 
floats,  this  means  that  the  node  is  held  fixed  in  the  vertical  direction 
but  free  in  the  lateral  directions.  All  three  components,  or  only  the 
vertical,  may  be  fixed  for  anchors.  Whenever  the  vertical  resultant  of 
all  of  the  element  tensions  connecting  to  the  point  exceeds  the  sum  of 
the  external  vertical  loads  at  that  point,  then  the  constraint  is  released. 
The  external  loads  are  assumed  to  include  the  distributed  loads  from  the 
elements  and  the  weight  or  buoyant  force  from  the  lumped  body. 

When  a  solution  step  allows  a  critical  node  to  move  past  the  limit 
position,  an  overshoot  procedure  is  activated.  On  incremental  solutions 
this  consists  of  a  step  division  in  which  the  portion  of  the  original 
step  which  will  satisfy  the  limit  and  the  remainder  of  the  step  is  done 
with  the  limit  restraint  imposed.  The  iterative  solution  imposes  the 
limit  by  moving  the  over-shooting  component  to  the  limit  state,  setting 
flags  to  prevent  immediate  release  of  the  constraint,  and  a  continuation 
of  the  iteration. 

4 . 7  Line  Payout  and  Reel -in 

The  treatment  of  line  payout  and  reel-in  presents  a  host  of  problems 
not  normally  encountered  In  structural  dynamics.  The  major  problem  is 
that  the  ma.ss  and  elastic  properties  change  with  time.  Not  only  does 
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this  produce  direct  changes  in  the  structure,  there  are  more  subtle 
changes  such  as  time  step  stability  characteristics.  For  example, 
during  reel-in  an  element  may  become  sufficiently  short  to  require 
reduction  of  the  time  step  to  obtain  convergence.  Other  problems  relat¬ 
ing  to  geometry  effects  will  be  pointed  out  as  this  discussion  proceeds. 

The  approach  taken  in  SEADYN  for  payout/reel-in  is  based  on  incre¬ 
mentally  changing  an  element's  unstretched  length.  Since  the  load 
carried  by  an  element  is  computed  from  element  strain  and  element  strain 
is  a  function  of  the  ratio  of  stretched  and  unstretched  lengths,  the 
equilbrium  state  is  unfluenced  by  this  length  change.  It  is  also  necessary 
to  adjust  the  mass  assigned  to  the  nodes  bounding  the  element.  Except 
for  some  procedural  bookkeeping  details  the  process  is  quite  simple. 

First,  it  is  assumed  the  payout/ reel-in  occurs  only  at  points  where  the 
motion  il;  defined.  It  may  be  fixed  or  given  specific  movements.  Second, 
the  incremental  change  in  element  length  is  computed  from  the  time  step 
and  the  average  payout  velocity  over  that  time.  The  node  where  the  pay¬ 
out  occurs  is  assigned  that  velocity  acting  in  the  instantaneous  direction 
of  the  element.  If  the  payout  point  is  a  moved  point,  the  velocity  of 
movement  is  added  to  the  point.  Mass  adjustments  are  then  made  to  the 
nodes  bounding  the  element  and  the  equilibrium  iterations  are  pursued. 

At  each  stage  of  the  iteration  the  velocity  components  at  the  payout 
point  are  recalculated  to  reflect  line  re-orientation. 

One  of  the  bookkeeping  problems  arises  when  payout  accumulates 
large  amounts  f  length  in  the  element,  or  reel-in  removes  a  major  portion 
of  the  clement.  In  cither  case,  a  process  is  initiated  which  modifies 
the  numoer  of  elements  and  nodes  in  the  system.  During  payout  the  modi¬ 
fication  involves  a  subdivision  of  the  element  and  the  addition  of  a  new 
node.  The  process  is  called  mitosis  for  obvious  reasons.  The  subdivision 
uses  linear  interpolation  to  locate  the  new  node  a  preselected  distance 
from  the  payout  point.  That  new  nodal  point  is  given  the  name  (number) 
of  the  payout  node  and  a  new  name  (number)  is  assigned  to  the  payout 
point  to  preserve  the  connectivity  specifications.  The  subdivided 
element  then  has  a  part  which  no  longer  grows  and  a  part  which  continues 
to  grow. 

The  reveise  process  is  still  calle<l  mitosis  even  though  that  is  a 
misnomer.  In  the  reel- in  process  the  element  which  shrinks  to  a  specified 
length  is  removed  from  the  system  by  adding  its  length  and  mass  to  the 
next  element  down  the  line  and  assigning  the  interfacing  node  to  the 
reel-in  point.  The  element  thus  removed  is  rendered  inactive  by  giving 
it  a  zero  unstretched  length.  The  node  removed  is  assigned  a  fixed 
status  to  avoid  solution  algorithm  problems. 

A  further  bookkeeping  problem  is  encountered  in  mitosis  when  the 
next  element  to  be  added  or  the  one  to  be  removed  is  not  of  the  same 
material  as  the  one  further  down  the  line.  This  multi-material  payout 
problem  is  dealt  with  by  an  approximation  as  outlined  below.  Consider 
a  payout  end  with  nodes  and  elements  are  inactive  and  awaiting  payout 
as  represented  in  Figure  4-1. 
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Figure  4-1,  Multi-material  Payout  Configuration  % 

At  the  beginning  cf  the  payout  sequence  the  active  payout  element  is  is 

entirely  composed  of  material  1.  The  new  length  added  to  that  element 
is  material  2  since  the  next  element  to  be  made  active  is  defined  to  be 
of  material  2.  Note  that  the  reserved  nodes  are  actually  given  the  saoie  _>v 

coordinates  with  fixed  conditions  on  all  components  and  the  elements  v 

connected  to  them  are  assigned  an  unstretched  length  of  zero.  > 


i 

I 

\ 


As  the  new  length  is  added  to  the  active  element,  two  lengths  are 
maintained  for  the  element.  The  first  is  the  length  at  the  beginning  of 
the  sequence  (i.e.,  the  length  of  material  1)  and  the  total  length  added 
(i.e.,  the  length  of  material  1).  Since  there  is  no  node  at  the  theoreti-  ^ 

cal  interface  the  active  element  is  treated  as  a  composite  element.  The  ;; 

composite  EA  and  effective  tension  are  calculated  as  follows:  :J 

1.  Use  current  nodal  positions  to  calculate  element  stretch,  i 

2.  Apportion  AL:  to  ALj ,  and  based  on  EA^ ,  and  EA2  from  the  ;; 

preceeding  time  step  or  iteration. 

3.  Calculate  T^,  T^,  EAj,  KA2»  material  properties  using  the 

new  sub-lengths  to  calculate  the  strains. 

4.  Compare  T  to  T„  and  repeat  steps  2  and  3  until  T  is  sufficiently 
close  to  T^. 

The  formulas  used  to  apportion  AL  to  materials  1  and  2  are  based  on  the 

small  strain  assumption  as  follovc:  ^1] 

Given: 


\  T,  EAj,  EA^ 
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Find: 


ot],  Qf,^,  ALj, 


Known : 


“^2  = 


t^, 


e  ~ 


AL 
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AE  *  T 

r  ^ 


small  strain  approximation 


Hence: 
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ihe  saiiit'  f(iuations  are  used  in  reol-ia.  The  reel-in  mitosis  will  assign 
the  length  of  the  removed  clement  (material  2)  to  the  next  element  down 
the  line  (material  1 )  by  a  A(.  equal  to  the  length  involved.  The  continued 
reel -in  will  then  remove  lengths  of  material  2  until  the  theoretical 
interface  point  is  reached. 

Yet  another  complication  is  encountered  when  lumped  bodies  are 
assigned  to  the  nodes  involved  in  payout/reci -  in .  When  a  lumped  body  is 
assigned  to  the  payout  node  the  theoretical  interface  is  again  employed. 
This  is  so  even  if  no  material  change  is  indicated.  Mass  is  then  shifted 
from  the  payout  point  to  the  first  active  node  as  line  is  paid  out. 

When  payout  mitosis  occurs  the  mass  is  reassigned  back  to  the  payout 
node  as  it  is  moved  out  to  the  new  active  position.  The  reverse  of  this 
process  is  employed  with  reel-in  of  lumped  bodies.  The  mass  redistribution 
formulas  are: 
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Where  represents  the  components  of  the  body  dragload. 

Finally,  another  bookkeeping  feature  has  been  found  to  be  necessary 
when  starting  a  payout  sequence.  Mitosis  is  defined  to  occur  when  the 
unstrctched  length  of  an  element  exceeds  a  user  defined  mitosis  length 
plus  a  reference  length.  The  reference  length  is  the  unstretched  length 
of  the  element  at  its  payout  initiation  or  the  mitosis  length.  The 
initial  length  is  taken  in  those  cases  where  there  is  raultimaterials 
involved  or  where  no  previous  mitosis  has  occurred  on  the  line  (i.e.,  at 
payout  startup).  The  reference  length  is  taken  to  be  the  mitosis  length 
in  all  other  situations.  This  multiple  choice  reference  length  is 
needed  to  avoid  premature  mitosis  when  starting  with  a  long  element,  and 
lo  retain  a  consistent  definition  of  the  theoretical  interface  on  multi¬ 
material  payout. 

A . 8  Drag  Aioplif icat ion  Due  to  Line  Strumming 

A  problem  of  considerable  importance  in  the  analysis  of  underwater 
cable  structures  is  the  estimation  of  the  forces  induced  by  interactions 
of  the  structure  with  the  moving  fluid.  The  assumptions  and  drag  models 
presented  earlier  neglected  an  important  but  ill-defined  phenomenon. 

When  the  cabie  cross  section  is  approximvitely  round  or  when  cable  fairing 
has  negligible  torsional  resistance,  the  steady-slate  value  of  the  lift 
force  is  negligible.  The  foregoing  developments  have  neglected  this 
force.  However,  under  appropriate  combinations  of  flow  conditions  and 
cable  size,  a  significant  oscillating  lift  force  can  be  pioduced  by 
vortex  shedding.  If  the  freqiu*ncy  of  this  oscillatory  lift  force  is 
sufficiently  close  to  the  frequency  of  a  structural  vibration  mode,  then 
significant  .structural  motions  are  possible.  Just  how  much  motion 
results  is  a  very  difficult  thing  to  predict. 

There  are  many  factors  unfluencing  the  structural  response  of 
flexible  systems  to  vortex  shedding.  Not  all  of  them  are  understood. 
Undoubtedly  such  things  as  structural  stiffness  and  orientation  (both 
local  and  global),  the  proximity  of  abrupt  changes  in  structural  features, 
the  coherence  of  vortex  shedding,  and  the  interaction  between  vortex 
formation  and  structure  movement  play  some  role  in  determining  the 
response.  A  comprehensive  theory  which  permits  complete  analysis  is  not 
yet  available.  Rather  than  wait  for  the  advent  of  a  complete  theory 
(and  in  true  engineering  tradition),  the  analyst  of  cable  structures 
introduces  an  approximation. 

In  many  situations  it  is  noticed  that  vortex  shedding  results  in 
relatively  small  amplitude  local  responses.  Rather  than  inducing  direct 
changes  in  global  structural  response,  the  cable  strumning  action  has 
indirect  effects  which  are  approximated  quite  well  as  an  apparent  increase 
in  cable  profile.  This  is  in  effect  an  increase  in  the  normal  drag  co¬ 
efficient.  Capitalizing  on  these  observations,  Skop,  Griffin  and  Ramberg 
|6)  of  NHL,  formulated  an  approach  to  computing  these  drag  amplifications. 
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Briefly  stated,  the  NRL  approach  presiimes  a  given  segment  of  cable 
which  is  a  candidate  for  strunming  will  struin  when  the  vortex  shedding 
frequency  is  in  proximity  of  one^gf  the  cable  natural  frequencies. 

Speci  f  ic.illy,  strutuaing  in  the  ii  cable  mode  is  assumed  to  occur  when 

u)  ^  u)  <  U)  . ,  4-10 

n  s  n+ 1 


where 


u)^,  are  two  adjacent  structural  natural  frequencies, 

u)^  is  the  vortex  shedding  frequency. 

This  presumes  that  the  natural  frequencies  are  closely  spaced  and 
tha  no  overlapping  in  responses  occurs.  When  lu  is  below  the  lowest 
natural  frequency,  it  is  assumed  that  no  strumming  occurs.  The  same 
assumption  is  made  when  lu  significantly  exceeds  the  highest  natural 
frequency. 

Once  the  critical  mode  is  identified  a  modal  scaling  factor  is 
computed  by 


,  -  ^  ^ 

n  “  i-L  ,  2  , 

where  (s)  represents  the  mode  shape. 

The  response  amplitude  along  the  string  is  estimated  by 
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The  value,  A  ,  is  a  dimensionless  modal  response  parameter  obtained 
from  a  least  squares  curve  tit  of  experimental  data  relating  it  to  a 
strumming  response  parameter.  S.. 


A  =  1.29/(1  +  0.43 

max  G 
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The  parameter,  S^,,  is  a  function  of  the  effective  damping  and  the  Vibra¬ 
tory  Reynolds  number.  The  latter  is  written 
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Where  u  is  the  kinematic  viscosity  of  the  fluid  and  d  is  the  effective 
cable  diameter.  The  strumming  response  parameter  is  given  by 


U-Vj 


“  2  n  K 

G  s 

where  S  is  the  Strotihal  number,  0.21,  and 

K  =  22.2/R  4-16 

s  n 

The  value  of  K  reflects  the  neglect  of  internal  cable  damping  and  the 
assumption  that  the  fluid  damping  coefficient  per  unit  cable  length  is 

=  4.5  n  p  u  R  4-17 

f  n 

where  p  is  the  fluid  density, 

The  strumming  response  estimate  of  Equation  (4-16)  is  then  used  to 
calculate  a  wake  response  parameter  as  follows 

W  =  (1  +  2  Y  ,  .)  (w/u»  )  4-18 

r  max/d"^  '  s 

Another  curve  fitting  of  experimental  data  leads  to  the  following  expression 
for  the  drag  coefficient  amplification 

C„(s)/C^  =  1.0  for  W  <1.0 

D  Do  r 

=  1.0  +  1.16  (W  -  for  W  ^  1.0  4-19 

r  r 

The  integrated  effect  of  the  local  drag  amplification  is  then  estimated 
by 

*^0 

The  procedure  requires  the  following  step.s: 

1.  Identify  a  specific  section  of  cable  which  is  a  candidate  for 
strumming. 

2.  Develop  an  approximate  local  model  for  the  cable  and  select  a 
representative  local  fluid  velocity. 

3.  Calculate  u) 

s 

=  2ns  v/d  I  Sin  4-21 
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whc  re 


V  =  roagni  tilde  of  local  fluid  velocity 
0^  =  velocity  heading 

0  =  representative  cable  heading 

A.  Calculate  natural  frequencies  and  mode  shapes  for  the  cable 
span. 

S.  Select  critical  mode. 

fc.  If  one  is  found,  calculate  drag  amplification  factor  using 
Kquations  (A-11)  through  (A-20). 

7.  Assign  this  drag  amplification  to  all  elements  of  the  string. 

This  procedure  is  repeated  each  time  it  is  determined  that  signifi¬ 
cant  changes  have  occurred  in  the  flow  or  structural  characteristics. 

A , 9  The  Initial  Configuration  Problem 

One  of  the  frustrating  features  of  cable  analysis  is  that  most 
cable  systems  do  not  have  rigidity  or  spatial  stability  unless  preloads 
are  imposed.  The  systems  are  usually  so  flexible  that  small  changes  in 
preloads  cause  large  shape  changes.  Noting  that  all  of  the  static  solu¬ 
tion  methods  use  a  stiffness  matrix,  one  is  faced  with  a  problem  of 
getting  a  realistic  estimate  of  the  stiffness  matrix  which  is  nonsingular. 
In  many  situations  it  is  necessary  to  have  very  accurate  estimates  of 
the  initial  configuration  before  any  of  the  solution  methods  will  work. 

Various  facets  of  this  problem  are  explored  in  Reference  1  and  7. 

Some  procedures  which  may  help  obtain  stable  starting  configurations  are 
discussed  in  the  User's  Manual. 

A  technique  that  has  proven  of  some  u.se  in  overcoming  ill-conditioning 
and  even  siugularilics  when  using  the  MNR  solution  is  the  use  of  numerical 
damping.  Felippa  |2ll  shows  that  noiisingular  ailjustment  to  the  estimator 
matrix,  jK)  can  be  generated  by  adding  a  matrix  of  the  form 

M  B  in  4-22 


where 


B  =  {R}  ’’iK]  {R}/{R}  4-23 

and  p  in  a  user  specified  numerical  damping  coefficient.  This  additional 
term  tends  to  "stiffen"  the  estimator  matrix  and  increase  the  chances 
for  convergence. 
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The  program  also  provides  for  numerical  damping  to  be  used  with  the 
liicr.'iiieiii  a  I  solutions.  This  leatuie  should  be  used  with  extreme  care 
since  it  alters  the  etpiations  of  eqiii  1  ihrium.  If  an  appropriate  value 
can  be  selected,  it  would  be  possible  to  compensate  for  some  of  the 
error  in  the  first  incremental  step.  There  is  no  rational  way  available 
to  estimate  how  much  damping  to  use  in  this  case. 

The  program  provides  a  quick  and  convenient  w.iy  of  getting  starting 
configurations  when  negatively  buoyant  lines  are  used.  If  a  line  between 
two  defined  points  is  a  catenary  which  reaches  the  bottom  with  a  horizontal 
tangent  somewhere  between  the  two  points,  then  nodes  can  be  generated 
along  that  line.  The  well  known  catenary  equations  are  used  and  one  or 
more  of  the  generated  nodes  are  assumed  to  be  on  the  bottom  if  the  lower 
defining  node  is  not  at  the  tangent  point.  The  element  preloads  for  the 
line  are  also  generated  and  can  be  assigned  to  the  line  elements. 

^.10  Component  Adequacy  Check 

The  SKADYN  program  provides  a  unique  feature  of  checking  the  capacity 
of  the  various  components  of  the  system  against  imposed  loads.  The  three 
types  of  checks  provided  are; 

1.  Anchors  -  the  loads  imposed  on  the  anchor  or  fixed  node  are 
summed  and  the  resultant  is  comp.'ired  with  the  anchor  holding 
power. 

2.  Buoys  -  the  resultant  of  the  load-s  in  the  lines  connecting  to 
the  buoy  is  checked  against  the  buoy  capacity. 

11.  Lines  -  individual  cable  elements  arc  checked  to  see  if  they 
are  loaded  beyond  their  capacities. 

The  componont  checks  tor  buoys  and  anchors  follow  the?  procedures  outlined 
in  NAVFAC  DM-26  and  discussed  in  Reference  22.  The  component  capacities 
can  be  input  or  obtained  from  the  inventory  tables  developed  for  the 
DKSMOOR  program  1221,  The  inventories  are  described  in  Appendix  D  of 
the  User's  Manual. 


All  of  lb','  adequacy  checks  rely  on  estimates  of  the  loads  at  the 
end  of  an  element.  T!ie  one  dimensional  simplex  eleoient  has  only  one 
v.ilue  of  tension  associated  with  it  regardle.ss  of  how  long  it  is,  what 
its  weight  is,  and/or  how  much  distributed  load  it  is  supporting.  The 
tension  associated  with  the  element  can  be  thought  of  as  being  at  the 
midpoint  of  the  element.  The  procedure  for  modeling  distributed  loads 
must  be  recalled  to  see  how  to  estimate  element  end  loads.  Equations 
(3-26)  and  (3-A3)  both  show  that  the  nodal  point  equivalent  loads  for 
distributed  load.s  arc  estimated  using  the  element  shape  functions. 

These  end  loads  represent  the  loads  applied  to  the  end  nodes  by  the 
element  as  it  supported  distributed  loads.  The  total  load  .it  each  end 
of  the  element  is  then  the  vector  sum  of  the  element  tension  and  the 
element  loads.  The  direction  of  the  resultant  load  at  each  end  gives  an 
estimate  of  the  direction  the  cable  lies  at  that  point.  In  symbolic 
form  the  total  loads  applied  to  the  nodes  by  the  element  are 
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Some  ambiguity  in  the  loads  applied  to  anchors  results  when  long  elements 
are  used  next  to  the  anchor.  If  the  line  actually  contacts  the  bottom 
before  tiie  anchor  connection  (i.e.,  between  the  two  nodes  which  define 
the  element),  the  SEADYN  Program  will  not  detect  this  in  the  solution. 
When  a  line  check  is  made  the  resultant  of  the  forces  at  the  anchors 
will  have  a  component  pulling  down  indicating  the  line  was  hanging  below 
the  .inchor.  This  is  due  to  the  model  not  being  able  to  sense  the  bottom 
contact  and  tr.')nsfer  the  weight  at  that  point.  Detailed  and  accurate 
modeling  of  the  line  interaction  with  the  bottom  requires  short  elements 
in  that  region.  Fortunately  the  leek  of  modeling  detail  at  the  bottom 
of  tlie  line  has  little  influence  on  the  response  at  the  top  of  the  line, 
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5.0 


SOLUTION  METHODS 


^  •  *  Static  Solutions 

The  SEADYN  program  offers  four  basic  solution  methods  for  static 
analyses.  The  various  features  of  each  of  them  are  discussed  at  length 
in  References  I  and  7  so  only  a  brief  description  will  be  given  here. 

The  first  method  is  a  .'.equence  of  linear  increments  (SLI  method). 

The  loads  are  divided  into  a  sequence  of  increments  and  the  basic  incre¬ 
mental  Equation  (3-12)  is  repeatedly  applied.  The  SLI  method  requires 
the  regeneration  of  the  incremental  stiffness  matrix  at  each  step  to 
reflect  the  changes  in  position  and  constitutive  relations.  It  has  the 
undesirable  feature  of  drifting  from  the  correct  solution  through  accumu¬ 
lating  errors  and  small  increments  are  required  for  accuracy. 

The  second  method  works  with  the  incremental  Equation  (3-13).  The 
first  step  is  identical  to  the  SLI  first  step,  but  on  each  succeeding 
step  the  force  residual  from  the  previous  step  is  fed  back  as  a  corrector. 
For  this  reason  it  is  called  the  residual  feedback  method  (RFB).  Although 
it  is  a  non-iterative  method  it  tends  to  be  self-correcting.  On  responses 
which  are  reasonably  well  behaved  (particularly  monotonic  responses)  the 
RFB  method  gives  accurate  results  with  significantly  fewer  steps  than 
does  the  SLI  method.  The  RFB  method  costs  somewhat  more  per  step  because 
of  the  residual  calculation.  The  recalculation  of  the  incremental 
stiffness  matrix  is  still  required  at  each  step. 

A  more  general  method  which  employs  iterations  to  solve  Equation  (3-11) 
is  known  as  the  modified  .Newton-Raphson  method  (MNR).  It  begins  with  an 
estimated  configuration  and  uses  an  estimate  of  the  tangent  stiffness 
matrix  to  obtain  successive  displacement  increments  which  hopefully  lead 
to  a  zero  force  residual.  Being  an  iterative  method  it  is  by  far  the 
most  accurate  of  the  three  methods.  This  accuracy  is  not  free,  however. 

In  the  first  place  the  method  is  not  unconditionally  convergent  and  the 
size  of  the  load  step  required  to  assure  convergence  is  not  easily 
determined.  In  some  cases  it  may  be  extremely  small.  In  ill-conditioned 
systems  it  may  not  be  possible  to  get  convergence  without  some  special 
auxiliary  procedures. 

The  general  form  of  the  MNR  method  is 

|K|  =  (R)“>  5-1 

where  the  superscript  i  refers  to  the  iteration  step.  The  (Kj  matrix  is 
referred  to  as  an  estimator  matrix.  If  (K]  is  the  incremental  matrix 
(K  1,  and  it  is  recalculated_at  each  iteration,  the  usual  Newton-Raphson 
procedure  is  obtained.  If  [KI  is  an  approximation  of  (ICj,}  and/or  it  is 
not  recalculated  at  each  iteration  one  has  a  modified  Newton-Raphson 
method.  When  (K}  is  not  changed  and  the  response  curve  is  monotonically 
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increasing  function  the  successive  estimates  will  usually  oscillate 
about  the  correct  solution.  If  the  step  size  is  small  enough  the  estimates 
will  tend  towards  the  solution  (converge).  If  the  step  size  is  too 
large  the  estimates  will  diverge. 

Once  oscillating  estimates  are  detected,  various  accelerating  pro¬ 
cedures  can  be  employed,  and  in  some  instances  they  will  work  even  when 
the  oscillations  are  divergent.  The  program  detects  oscillations  by 
monitoring  the  degree  of  freedom^ylth  the  largest  initial  response 
(i.e.,  largest  component  of  {Aq}^  .  Oscillation  is  signalled  by  repeated 

sign  changes  of  the  critical  {Aq}  component.  The  sisiplest  acceleration 
scheme  averages  the  two  alternating  estimates  using  the  sizes  of  the 
critical  components  to  weight  the  average.  An  optional  procedure  uses  a 
one-dimensional  search  to  seek  a  close  estimate  of  {Aq}  which  crosses 
the  correct  solution.  This  search  is  initiated  when  the  ith^^teration 
signals  a  large  oscillation.  The  search  begins^gt  the  (i-1)  position 
and  takes  small  steps  in  the  direction  of  the  i^  increm»'t.  A  new  step 
is  tried  at  each  position  until  the  new  increoient  reverses  direction. 

At  this  point  the  last  two  positions  are  averaged  to  get  a  new  starting 
estimate  for  the  iteration.  The  size  of  the  step  used  in  this  1-D 
s^^rch  is  controlled  by  input.  The  input  factor  is  the  fraction  of  the 
i''  step  which  is  taken  as  the  first  step  of  the  search. 

Various  options  are  provided  to  measure  convergence  of  the  itera¬ 
tions.  Reference  1  should  be  consulted  for  discussions  of  the  convergence 
criteria.  The  options  are  listed  in  Sections  7.1.14  of  the  User’s  Manual. 

The  fourth  solution  method  is  called  a  viscous  relaxation  method. 

It  is  similar  in  '".'r.y  respects  to  the  damped  Newton-Raphson  method.  The 
major  distinction  xs  that  the  equations  are  cast  into  a  transient  dynamic 
format  where  a  pseudo-time  is  used  to  control  load  steps  and  iterations. 

A  pseudo-velocity  vector  is  used  in  conjunction  with  the  residual  vector 
to  make  modifications  in  tne  damping  and/or  time  steps.  This  allows  the 
solution  to  be  adapted  to  the  character  of  the  response.  As  damping 
becomes  negligible,  the  method  becomes  the  Newton-Raphson  method.  Details 
of  the  derivation  are  giver  in  Reference  7.  The  solution  form  is 

(sit  =  r"^')  -  ‘f’f  5-2 

*  (¥) 

Approximating  at  the  geometric  positions  of  linearizes  the 

method.  It  is  applied  the  same  way  as  the  RFB  method  except  that  the 
time  step.  At,  and  the  damping  matrix  are  adjusted  based  on  the  process 
of  the  solution.  The  response  data  for  the  step  are  cosiputed  from 

+  At  ^'{q}  ♦  oAt  -  *'{q}  5-3 
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Th«>  (lamping  matrix  is  computed  from 

(diagonal ) 
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where 

is  a  modifiable  damping  parameter 
Cj  is  a  response  dependent  term 
K  is  a  representative  stiffness 

The  value  of  K  is  taken  to  be  the  average  extensional  stiffness  (AE/L) 
or  Equation  (4-23).  The  C.  term  is  included  for  special  problems  associated 
with  the  angular  responses^of  rigid  bodies. 


The  pattern  of  velocity  changes  and  the  residual  norm  are  used  to 
adjust  At  and  ^  in  controlling  the  performance  of  the  algorithm. 

The  form  using  Equation  (5-2)  is  referred  to  as  the^yW  metho^.  A 
form  which  avoids  the  residual  calculation  by  replacing  {f}  -  {g} 

by  {Af}  is  referred  to  as  the  VRS  method.  Its  performance  is  inferior 
to  the  VRR  method. 


5 , 2  Time  Domain  Dynamic  Solutions 

The  SEADYN  program  provides  for  four  different  time  domain  solution 
method.s.  Each  of  them  utilize  the  following  set  of  forward  difference 
equations  to  develop  the  solutions: 


+  At  ‘'{q}  +  ^  *-{q}  +  a  At  -  *"{4})  5-7 

2 

+  PAt^  (^^^‘^{q*}  -  ^{q]) 

=  ^-Iq}  +  At  ^{ql  +  TAt  -  ^{q]) 

where  At  is  the  time  step  and  a,  p,  and  T  are  integration  parameters. 

The  usual  Newmark  formulas  are  obtained  with  a  =  0. 
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When  Equations  (5-7)  are  augmented  with  a  motion  equation  for  the 
system,  one  has  three  simultaneous  dif ferential^^gyations  with  the  nodal 
displacements,  velocities  and  accelerations  at  C  as  unknowns. 
Specific  forms  of  the  solution  routines  for  this  set  of  equations  are 
developed  in  Reference  1.  Only  the  final  forms  will  be  presented  here. 

Direct  Iteration  Method  (DIM) 

This  method  uses  Equations  (5-7),  (3-8),  and  (3-46)  to  obtain  an 
iterative  formulation.  The  form  is 


(„*1)  ^  (1,),  5-8 

l.At,^|(r,.l)  ^  ,  TAt  A(""|A^j '"*’>) 

The  iteration  can  be  started  with  fJquations  (5-7)  with  a  =  p  =  T  =  0  or 
with  tiie  residual  feedback  solution  to  be  described  below.  The  DIM 
method  has  been  shown  to  be  an  accurate  and  cost  effect  approach.  Since 
it  is  an  iterative  method  which  retains  all  of  the  nonlinear  terms  it 
readily  deals  with  all  of  the  important  nonlinear  effects  such  as  position 
depesident  loads,  constraints,  slack  segments,  and  changing  geometry. 

Of  particular  interest  is  the  ease  with  which  problems  with  defined 
motion  arc  dealt  with.  When  problems  require  the  imposition  of  a  known 
motion  (e.g.,  cable  towing  problems)  the  moved  node  is  simply  held  at 
the  required  dynamic  state  during  the  iterations. 

of  Linear  increments  (S1.1) 

This  method  starts  from  Equations  (5-7)  and  uses  the  incremental 
form  (3-9).  An  iterative  solution  is  also  employed,  but  it  has  been 
shown  that  this  method  costs  essentially  the  same  .is  the  DI  method  and 
IS  !im<  ii  less  accurate  (11.  It,  loo,  is  capable  of  solving  the  moving 
bi.'aniiary  problem  but  an  augmented  global  stiffness  ma’.rix  involving  the 
moved  degrees  of  freedom  is  required.  This  increases  the  storage  require¬ 
ments,  reduces  the  speed  and  increases  the  xiumerical  error  potential. 

Ke.^idua  1  feedback  Method  (RfB) 

Eqiiations  (5-7)  are  inverted  and  substituted  into  the  incremental 
form  (3-10).  The  result  is 


;  ;f 


1  - 


5-9 
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where 


=  (K  )  +  .-i—  {M(''{q}))  5-10 

=  -  Mg}  ^  lM(Mq}))  5-ii 

-  (1  -  ^)  Mq}) 

rquat-ion  (5-9)  is  a  linear  algehraic  equation  which  can  be  solved  by  the 
usual  manner.  No  iteration  is  required  and  the  method  has  a  self-coriecting 
feature  .similar  to  the  static  KFB  method.  This  method  follows  more 
closely  tlie  traditional  Newniark  implicit  integration  form  [23].  It  has 
the  unconditional  stability  features  that  have  made  the  method  so  popular. 

No  provision  is  made  for  solving  the  moving  boundary  and  payout  problem.?, 
however.  The  fact  that  the  RFB  method  requires  the  formation  of  a 
global  stiffness  matrix,  the  evaluation  of  .i  re.sidual,  and  the  solution 
of  s iinul  t .sneous  equations  reduces  its  cost  effectiveness.  This  is 
■  sually  compensated  for  by  using  larger  time  ste^/s. 

.qodi_fi<.'d  Nevdon-Raphson  Forin^ 


IK 


eff ' 


rtns  IS  an  iterative  method  based  on  a  procedure  similar  to  the  MNK 
.static  solution.  The  form  is 


IK^.l 
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Tlif  Iteration  is  .started  will: 
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where  {Aq}  is  obtained  from  an  RFB  solution.  As  an  alternative 
to  Equations  (5-13)  the  iterations  can  be  started  from  Equations  (5-7) 
with  a  =  p  =  T  =  0. 


The  accuracy  of  this  method  has  been  demonstrated  but  its  cost  per 
time  step  does  not  appear  competitive  with  the  DI  method  [1].  Moving 
boundary  and  payout  solutions  have  not  been  implemented  with  this  solution. 


It  would  appear  from  the  studies  done  in  Reffcr..'nce  1  that  the  best 
general  purpose  solution  method  of  the  four  is  the  .'Jl  method.  On  some 
problems  the  method  may  encounter  some  difficulties,  however.  The 
Iterations  are  convergent  only  for  step  sizes  of  the  o.-der  it  takes  an 
axial  wave  to  traverse  the  shortest,  stiffcst  element.  When  the  imposed 
forces  arc  slowly  varying  the  routine  will  attempt  to  increase  the  step 
size  beyond  stable  behavior  limits,  because  the  motions  are  small  in 
small  time  steps  and  the  iterations  converge  rapidly.  Controls  are 
provided  to  the  user  to  prevent  step  size  increases  in  this  situation. 

In  some  situations  the  RFB  solution  may  be  found  very  cost  effective. 
Generally  this  will  be  nearly  linear  systems,  which  are  slowly  varying. 

The  user  is  cautioned  against  using  very  large  steps  with  the  RI'B  method 
since  gross  errors  may  result  (1). 

it  should  be  noted  that  the  diagonal  form  of  the  mass  matrix,  [M^)> 
is  used  in  each  of  the  above  methods.  The  advantages  in  storage  and 
solution  effort  are  obvious.  It  has  been  shown  that  with  ptOpt;r  selec¬ 
tion  of  the  integration  parameters,  the  errors  introduced  by  the  lumped 
mass  approximation  tend  to  compensate  for  those  induced  by  the  solution 
algorithm  to  give  more  accurate  results.  Experience  has  shown  that 
0=0,  P  “  1/12  and  T  -  1/2  is  good  choice  of  parameters.  Preliminary 
observation  (1)  indicate  that  a  small  negative  value  for  Of  may  be  benefi¬ 
cial  ill  controlling  amplitude  atteniiati on.s .  Values  of  T  >  1/2  introduce 
numerical  damping  of  the  higher  frequency  components  while  T  <  1/2 
introduces  negative  damping. 


5 . 3  Frequency  Domain  Solution 

5.3.1  Solutions  for  Regular  Waves 

The  response  amplitudes,  (Q) ,  for  a  given  frequency  are  obtained  by 
solving  the  linear  simultaneous  algebraic  equations  represented  by 
(.3-16).  The  coefficient  matrices  for  mass  and  damping  must  be  recalculated 
for  eacti  frequency  since  the  linearization  procedure  leaves  them  dependent 
on  the  wave  frequency. 
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The  damping  terras  present  some  further  difficulty.  In  addition  to 
being  frequency  dependent,  the  lineari/.ed  viscous  terras  are  dependent  on 
the  amplitude  of  the  response.  The  ship's  roll-damping  depends  on  the 
roll  angle,  the  buoy  rotational  damping  depends  on  the  rotation  angle, 
and  the  cable  and  lumped  body  damping  depend  on  the  lateral  displacement 
amplitudes.  Thus,  it  is  seen  that  the  incremental  equations  a  'e  not 
strictly  linear. 

An  approximation  procedure  is  introduced  to  deal  with  this  problem. 

This  involves  iterative  solutions  of  the  Equation  (3-l6)  for  each  frequency. 
The  first  solution  at  a  given  frequency  is  calculated  assuming  a  ship's 
roll  angle.  The  roll  angle  obtained  from  the  solution  is  then  used 
along  with  the  other  pertinent  response  amplitudes  to  recalculate  the 
damping  terms  and  obtain  another  solution.  This  procedure  is  repeated 
until  two  successive  estimates  of  ship's  roll  are  within  1“  of  each 
other.  It  is  assumed  that  buoy  ;'.nd  cable  damping  are  less  important 
than  ship's  roll  and  are  thus  converged  when  the  roll  has  converged. 

The  response  is  then  dependent  not  only  on  the  wave  frequency,  but  also 
on  the  wave  amplitude.  This  means  that  it  is  not  appropriate  to  assume 
a  unit  amplitude  for  a  given  wave  frequency  with  the  intent  of  obtaining 
a  Response  Amplitude  Operator  (RAO)  for  that  wave.  It  is  necessary  to 
have  the  correct  wave  amplitude  at  each  frequency.  Therefore,  the  sea 
,spect.'um  must  be  used  in  the  calculation  of  the  regular  wave  responses. 

Once  the  steady  state  response  for  a  given  wave  frequency  and  amplitude 
IS  oblai'ied,  the  RAO  is  estimated  by  dividing  by  the  wave  amplitude. 

The  program  allows  the  mass  matrix  to  be  formed  either  with  the 
lujnped  or  consistent  form  of  the  clement  mass.  The  lumped  form  is  used 
for  the  fluid  added  mass  terms. 

Provisions  for  iut«?rnal  damping  effects  are  provided  in  the  propor- 
lion.il  damping  form.  Thus 

IC]  =  a  (M]  +  p  (K,j.]  5-14 

where  a  and  p  are  proportionally  constant.s.  The  damping  from  fluid  drag 
effects  on  the  cable  elements  most  be  linearized  before  it  can  be  used 
in  Equation  (3-16).  This  linearization  is  described  in  Appendix  I. 

5.3.2  The  Steady-State  Wave- Induced  Drift  Forres 

Whenever  waves  encounter  a  floating  body  there  arises  a  set  of 
forces  which  tend  to  move  the  body  in  the  direction  the  wave  is  traveling. 
These  forces  are  often  neglected  since  they  are  usually  small  in  magni¬ 
tude.  These  are  the  so-called  second  order  wave-induced  drift  forces. 

They  are  generally  slowly  varying  compared  to  the  frequency  of  the  inci¬ 
dent  wave;  however,  they  have  an  average  or  steady-state  component  which 
may  be  .significant  enough  to  cause  an  adju.stment  of  the  static  position 
of  the  ship.  These  forces  are  directionally  dependent  and  are  sensitive 
to  the  amplitude  of  the  ship’s  response  to  the  wave.  The  DTNSRDC  Ship's 
Motion  File  provides  a  table  of  coefficients  which  can  be  used  to  estimate 
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the  steady-state  drift  forces  after  the  ship’s  dynamic  response  is 
obtained.  It  should  be  noted  that  these  forces  do  not  estimate  the 
dynamic  effects,  which  are  at  a  lower  frequency  than  the  incident  wave. 

The  specific  form  for  the  drift  forces  is  given  in  Appendix  A  of 
the  User's  Manual  in  the  discussion  of  the  Ship  Motion  File.  The  result 
is  a  set  of  forces  for  the  lateral  and  longitudinal  directions  and  a  yaw 
moment  acting  at  the  ships  reference  point  for  each  wave  frequency.  It 
is  assumed  that  the  drift  forces  are  cumulative  for  the  various  wave 
represented  in  a  wave  spectrum.  Therefore,  the  drift  forces  for  each  of 
the  regular  waves  are  accumulated.  It  is  felt  that  this  is  a  reasonable 
approximation  since  the  wave  amplitude  indicated  by  the  wave  spectrum  is 
used  in  response  calculation. 

These  accumulated  drift  forces  c^n  then  be  used  as  an  additional 
static  loading  to  adjust  the  static  reference  state.  If  it  is  felt  that 
this  adjustment  will  affect  the  regular  wave  solutions  significantly, 
the  user  can  request  iterations  on  the  regular  wave  solutions  and  the 
drift  force  adjustments  to  the  static  reference  until  sufficiently  small 
changes  are  found.  Either  the  MNR,  VR,  or  RFB  solutions  can  be  used  to 
adjust  the  static  reference. 

5.3.3  Solution  Procedure  for  Random  Seas 

The  superposition  of  the  regular  wave  responses  to  represent  the 
response  to  random  seas  follows  the  well  established  methods  from  the 
theory  of  random  vibrations  [2Al.  The  sea  is  assumed  to  be  uni-directional 
(long-crested)  and  is  described  in  terms  of  a  generalized  spectral 
energy  density  function  having  the  following  form: 


S(u)) 


.  ,  5  -B/w 
A/u)  e 


5-15 


Typical  values  for  the  parameters  A  and  B  are  given  in  Section  7.3.2  of 
the  User's  Manual. 


The  frequencies  to  be  included  in  the  set  of  regular  wave  calculations 
<nre  determined  by  specifying  a  frequency  increment,  iiw,  a  lower  bound 
and  upper  bound,  Regular  wave  responses  are  than  calculated  for 

u),  =  u)  .  +  ^  Aio  +  (i-l)Aw 

1  min  2 


At  each  frequency  the  incident  wave  height  is  determined  from  the  sea 
spectrum  by  the  following: 

h(ut.)  =  (|s(ui.)Aiu)^^^  5-16 

(Note  that  the  wave  height  is  twice  the  wave  amplitude.)  This  wave 
height  is  then  used  to  converge  on  a  steady-state  dynamic  response  using 
the  methods  described  in  the  previous  section. 
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Let  H(u),)  represent  the  response  of  one  of  the  quantities  (nodal 
displacement  component  or  element  tension).  The  response  spectral 
density  of  this  quantity  is  then  given  by 

S^(u))  =  HH*  S(u»)  5-17 

The  mean  square  of  the  response  is 

CB  oa 

F.(x^]  ~  / HH*  S(iu)du»  =  /  S  (uj)du)  5-18 

0  0  ^ 

The  integral  is  evaluated  numerically  using  the  points  obtained  from 
each  of  the  regular  wave  solutions.  Thus 

9  N 

Elx^l  ~  ^  S^0-^)Au)  5-19 

i=l 

where  N  is  the  number  of  regular  wave  components  used. 

When  the  response  quantities  represent  the  dynamic  excursions 
relative  to  the  static  reference  state,  then  their  expected  values 
(i.e.,  their  means)  are  zero.  The  magnitude  of  the  response  is  then 
treated  as  a  static  part  plus  a  dynamic  part.  The  amplitude  of  the  dy¬ 
namic  part  is  assumed  to  follow  a  Rayleigh  distribution  with  a  mean-square 
value  which  is  a  function  of  the  area  under  the  response  spectrum  curve. 
Assuming  the  sea  spectrum  used  is  based  on  double  the  square  of  the  wave 
height,  the  mean-square  of  the  response  amplitude  is  the  value  obtained 
(.5-19)  divided  by  eight.  It  is  possible,  then,  to  make  statistical 
estimates  of  the  maximum  response  by  making  statistical  estimates  of  the 
dynamic  part  and  adding  them  to  the  static  values  obtained  in  the  updated 
reference  configurations. 

5  .A  Njr  t  i^rjj  1  F  r  e  queue  i_e  ^  ^iliL 

The  incremental  Equations  (3-9)  without  damping  or  external  forces 
has  the  form 


IM)  {u)  +  [K.j.j  {u}  ^  0  5-20 

If  it  is  assumed  that  simple  harmonic  motion  occurs  then 

(-lu^  (M]  +  (K.J,]  {u}  sin  uit  =  0  5-21 

since  {u}  sin  u)t  =  0  only  in  the  trivial  case  of  {u}  =  0,  this  requires 
that 


det  (jK,j.|  -  (Ml)  =  0 


5-22 
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Solutions  to  this  equation  lead  to  the  natural  frequencies  and  mode 
shapes  of  the  small  displacement  free  oscillations  about  the  configurations 
used  in  expressing  [M]  and  [IC] .  An  iterative  solution  is  used  which 
obtains  all  of  the  natural  frequencies  and  mode  shapes  for  the  system. 

The  primary  assumption  made  in  getting  the  frequencies  and  mode 
shapes  is  that  the  equations  are  linear  and  no  significant  damping 
exists.  Both  of  these  assumptions  may  be  violated  to  some  extent  for 
underwater  cable  structures.  The  option  is  provided  since  the  frequencies 
and  mode  shapes  may  still  be  indicative  of  the  structural  behavior. 

5 . 5  Time  Sequenced  Static  Solutions 

Often  the  transient  response  problem  involving  boundary  movement 
and/or  line  payout  can  require  considerable  computer  resources  to  evaluate. 
In  situations  where  imposed  velocities  and/ur  payout  rates  are  relatively 
slow,  approximate  solutions  can  be  obtained  by  a  sequence  of  static 
"snap  shots."  SEADYN  has  implemented  an  option  to  drive  the  static 
solutions  sequentially  through  boundary  motion  and  payout  conditions. 

This  option  is  called  the  Time  Sequenced  Static  Solutions  (TSSS).  Time 
is  simply  a  parameter  in  this  case  which  determines  where  the  moving 
boundary  is  and  how  long  the  unstretched  payout  elements  are.  Only  the 
iterative  static  solutions  (MNR,  VRR,  VRS)  are  valid  in  this  option. 

The  previous  state  is  used  as  a  starting  guess  for  each  iteratively 
solved  "snap  shot."  Time  steps  can  be  as  large  as  the  starting  guess 
will  allow;  i.e.,  as  long  as  the  previous  state  remains  a  reasonable 
guess  for  the  next  state. 
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